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Abstract 

Collective interactions in the preequilibrium quark matter and 
hadronic resonance gas stage of ultrarelativistic nucleus-nucleus col- 
lisions are studied in the framework of the the transport theoretical 
approach RQMD. The paper reviews string fusion into color ropes 
and hadronic rescattering which serve as models for these interactions. 
Hadron production in central Pb(160AGeV) on Pb collisions has been 
calculated. The changes of the final flavor composition are more pro- 
nounced than in previous RQMD studies of light ion induced reactions 
at 200AGeV. The ratio of created quark pairs ss/ {uu+dd) is enhanced 
by a factor of 2.4 in comparison to pp results. Color rope formation 
increases the initially produced antibaryons to 3 times the value in the 
'NN mode', but only one quarter of the produced antibaryons survives 
because of subsequent strong absorption. The differences in the final 
particle composition for Pb on Pb collisions compared to S induced 
reactions are attributed to the hadronic resonance gas stage which is 
baryon-richer and lasts longer. 
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1 Introduction 



The ground state of quantum chromodynamics shows pecuhar properties, 
confinement of colored degrees of freedom and spontaneous breaking of chi- 
ral symmetry. It is expected that at high temeratures and large baryon 
densities chiral symmetry gets restored and quarks are no longer confined. 
The central goal of modern heavy ion physics is to explore these properties 
of strongly interacting matter |1|, 0. Several experiments on nuclear tar- 
gets are performed nowadays using the Au(11.6AGeV) beam at BNL or the 
Pb(160AGeV) beam at CERN. The first series of experiments with light ion 
beams (p,0,Si,S) has just been finished in these laboratories. 

The work on which this paper is based has started some years ago [^]. 
The goal was set to develop a complete transport theoretical scenario of 
nucleus-nucleus reactions, from the initial state of two nuclei before overlap 
to the final state after the strong interactions have ceased (freeze-out). The 
developed model has been dubbed relativistic quantum molecular dynamics 
(RQMD). RQMD is a semi-classical microscopic approach which combines 
classical propagation with stochastic interactions. Strings and resonances can 
be excited in elementary collisions. Their fragmentation and decay leads to 
production of particles. Overlapping strings do not fragment independently 
from each other but form 'ropes', chromoelectric fiux-tubes whose sources 
are charge states in higher dimensional representations of color SU^ @]. The 
nature of the active degrees of freedom in RQMD depends on the relevant 
length and time scales of the processes considered. In low energy collisions 
(around 1 AGeV) RQMD reduces to solving transport equations for a system 
of nucleons, other hadrons and eventually resonances interacting in binary 
collisions and via mean fields (similarly to BUU [Q, QMD 0, and so on). 
At large beam energies (> 10 AGeV) the description of a projectile hadron 
interacting in a medium (in the simplest case a cold nucleus) as a sequence of 
separated hadron or resonance collisions breaks down 0. A multiple collision 
series can be formulated on the subhadronic (quark) level. Following the 
paths of the ingoing constituent quarks, a Glauber-type multiple collision 
series is generated in RQMD, with cross sections taken from the additive 
quark model . The secondaries which emerge from the fragmenting strings, 
ropes and resonances may interact with each other and the original ingoing 
hadrons (rescattering and mean- field interaction). 

The RQMD model has been successfully applied to nuclear reactions at 
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ultrarelativistic energies (see e.g. 0-|T^). The use of transport models like 
RQMD is not restricted to the study of the generated final state. It is of 
equal importance that such models can be used to study the influence of 
various types of interactions and medium effects on final state observables 
and to achieve a better insight into the transient stages of heavy ion colli- 
sions. Using the information from calculated RQMD events these questions 
have been addressed already in the literature, at 200AGeV mostly for S on 
A collisions |11]-|14]. Here and in a follow-up paper |15] I wish to expand on 
those results and study the amount of collectivity for central Pb(160AGeV) 
on Pb collisions. This is the heaviest projectile-target combination for which 
experiments are currently undertaken at CERN-SPS. Any kind of collec- 
tivity - physics of dense matter beyond mere superposition of independent 
nucleon-nucleon interactions - is expected to be strongest in a system with 
minimized surface-over-volume ratio. The paper presented here will focus 
on the influence of some important building blocks of RQMD - color ropes 
and hadronic rescattering - on the particle chemistry and distributions in 
phasespace. In particular, I shall consider the yields and momentum distri- 
butions of strange hadrons and antibaryons. Their production is dynamically 
suppressed in elementary hadronic interactions which makes them a useful 



probe for collectivity in the transient stage of AA collisions |]I6| . The follow- 
up paper mainly discusses the dynamical evolution of the created source in 
space-time until freeze-out. 

The achieved degree of collectivity in the first - the prehadronic or quark 
matter - stage of nucleus-nucleus collisions is a central topic of heavy ion 
research. Soft multiparticle production is commonly described by fragmen- 



tation of excited color strings (see for reviews on this subject |T^, ^ |2T|). It 
is expected that independent string fragmentation, which means no collec- 
tivity at all, breaks down in central nucleus-nucleus collisions, because the 
string density gets too large. Interactions of overlapping strings is modeled 
as color rope formation in RQMD. Ropes can be viewed as a model of locally 
deconfined quark matter, which is dominated by longitudinal excitations and 
therefore of relevance for the preequilibrium stage in nucleus-nucleus colli- 
sions. The evolution in a rope is entirely governed by nonperturbative dy- 
namics (a distance-independent rope tension). Other 'second generation' 
transport approaches also go beyond independent string fragmentation, for 
instance string fusion in QGSM [021, the Spanish version of the dual-parton 



model [23 1, or quark matter droplet formation in VENUS [24|. A common 
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result of these approaches is that strangeness and baryon pair production 
which is dynamically suppressed in elementary reactions becomes more fa- 
vorable. On the other side, arguments have been put forward to include 
diquarks and an unsuppressed strange quark component into the nucleon sea 
'from the beginning' which can be transformed into real particles by multiple 
collision effect s[^ . 



Collectivity may also emerge in the hadronic stage of nuclear reactions, 



resonance matter formation - as studied in the models RQMD [^, and 
ARC - or mean fields [^, In recent years the interactions of reso- 
nances in a dense medium have found a lot of attention. For instance, res- 
onances may act as 'energy storage' in multi-step collision processes which 
are of importance for heavy particle production |31|, Strictly spoken, 
there is no scattering (S) matrix for resonances, because they are unsta- 
ble. However, introducing them as quasi-particles makes sense under some 
limiting conditions (F/M <^ 1). Note that resonance interactions can get 
important only in a system at sufficient density. The collision frequency has 
to be larger than the inverse of resonance lifetime (typically 1-2 fm/c). 

The RQMD approach has evolved over the last years, because various 
interaction pieces have been put in step by step. After the important role 
of resonances had been recognized, a model was constructed for annihila- 
tion of mesons and baryons, including resonances, into nonstrange baryon 
p6| and into a 'string continuum' [^]. Independent string frag- 



resonances 



mentation was replaced by color rope formation in case of overlapping strings 



np. Low energy hadronic interactions - s channel resonance formation, t 
channel meson and Reggeon exchange - have been modeled in the meson- 
meson |60] and the meson-baryon sector with net strangeness [14|. In this 



paper we are going to review these interactions in some detail. Some recent 
work to further develop the RQMD model for applications up to collider en- 
ergies (multi-string excitations according to the AGK cutting rules [Q) will 
be described elsewhere |^5|]. These interactions have been implemented into 
the computer code RQMD (most recent version 2.1) which was developed by 
the author. 



^ For the topics of interest here multi-string excitations in NN cohisions are of minor 
importance. On the one side, they tend to enhance the string densities. However, with 
respect to rope formation most of this effect is canceled at the comparably small CERN- 
SPS energy, because the average string length decreases in comparison to the case with 
2-string excitations only. 
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2 Strings and ropes in hh and AA collisions 



The so-called recombination approaches like RQMD assume that in elemen- 
tary hh interactions the hadron wave function is decomposed into spectator 
and interacting quarks. The quark spectators neutralize their color while 
keeping their gluon cloud coherently bound (constituent quark picture). See 
Ref. for a more detailed description. The interactions are dominantly 
initiated by slow quanta. The interacting constituent quark looses part of its 
gluon cloud, because the interacting components get out of phase with the 
spectator remnants in the cloud. These gluons form - together with their 
partners from the target - the source for secondaries. The parton language 
is used here rather loosely, to give a rough sketch of the process. There will 
be no attempt here to describe this process quantitatively by interactions 
of QCD quanta. It will be simply assumed that these interactions result 
in the longitudinal excitation of color strings. The concepts of strings in 
strong interactions actually predates QCD. There are many recent attempts 
to understand the properties of strings or equivalently chromoelectric flux- 
tubes directly from QCD, e.g. by invoking the dual Higgs mechanism for 
QCD-monopole condensation fl^ . 



2.1 Strings in RQMD 

In inelastic hh collisions the fraction of the projectile lightcone momentum 
P"*" = E + pi {pi > assumed) which goes into target string excitation is 
determined from 

, (1) 

with the lower limit set by the target hadron's (small) plus-component before 
collision. The target momentum (negative lightcone component) which goes 
into projectile string excitation is determined in the same way. This string 
excitation law is the same as originally suggested by the Lund group and 



realized in the Monte Carlo code FRITIOF |T^. The invariant masses of 
the target and projectile excitations should not fall below some minimum to 
allow for particle production. The minimum value is set to m+0.3 GeV/c^, 
with m usually taken as the ingoing hadron mass. Since the ingoing pseu- 
doscalar mesons which are Goldstone bosons have exceptionally light masses, 
m is set in these cases to 0.7 (0.85) GeV for mesons without (with) strange 
valence quarks. The excited states decay either stringlike or - for excited 
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masses below some 'smeared-out' cut-off value (2 GeV/c^ for nucleons) - as 
resonances. 

The string topology generated in RQMD for hh interactions and displayed 
in fig. |l| is determined from the basic assumptions that all ingoing valence 
quarks keep on moving into the same direction and no net color flows between 
target and projectile. A sea quark pair gets polarized in course of the inter- 
action resulting in spatial separation of sea quark charge and anticharge (see 
fig. m). The interacting constituent quark combines with the sea anticharge 
into a colorless state. In the semi-classical picture the sea-quark compan- 
ion starts to move backward due to the momentum flow from the target. 
Confinement forces the backward moving quark to pull out a string. The 
spectator quarks neutralize color in the same way as the interacting quark 
at the other end of the string. Of course, if energy is lacking the backward 
and forward sea-quark pair might be the same. In this case there will be no 
string, but decay of the created hadron excitation into two hadrons. In the 
terminology employed here a 'color string' (restricted to be a longitudinal 
excitation) is spanned between sea quarks which are identified with partons. 
Such a configuration is assumed to fragment in the same way into hadrons 
as a qq configuration produced in e+e~ annihilation. In contrast, the ingoing 
valence quarks are assumed to be 'dressed' (constituent quarks) and better 
prepared for hadronization. An ingoing spectator quark fragments differ- 
ently (harder) than a partonic quark with equal momentum giving rise to 
the so-called leading particle effect. In the model quark spectators are able 
to transfer all their primordial momentum to the hadron(s) into which they 
fragment. The corresponding constituent quark structure functions are spec- 



ified in [|T3|, |3^. The interacting quark keeps only a fraction of its original 
momentum, in contrast to the spectator quark(s). The momentum which 
it keeps is dubbed in fig. |l]. The difference to the original momentum 
PyoYo is used up as the forward momentum of the excited string, i.e. for par- 
ticle creation, is determined stochastically using the string fragmentation 
function f{z) [jl9| as the probability measure. The backward momentum of 
the string p^^ transferred from the target is determined from the equivalent 
of eq. (|ip. The decay properties of a color string are completely determined 
by its light cone momentum and the flavor at its end. The spectator mo- 
mentum pI' does not enter at all into the calculation of the string evolution 
or its decay products, (see fig. 0). 

The decay of elementary color strings into hadrons is calculated employ- 
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ing the concept of 'left-right symmetric string fragmentation' developed by 



the Lund group ||T9[, with default JETSET 6.2 default parameters for f{z) 
[ p^ . The formation points of hadrons from string decay are calculated as 
the average of the two break points from which the quark constituents are 
emerging. This is schematically displayed in fig. 1 from which one can also 
read off the prescription for the formation points of hadrons containing one of 
the original valence quarks. In particular, the formation length of the lead- 
ing hadrons coincides with the concept of so-called 'constituent formation 
length' [^1 in which the formation length shrinks to in the extreme limit 
z—>-l. The standard high energy 2-string scenario is modified in three cases, 
generation of additional strings from sea excitations which can be related to 
multi-Pomeron-exchange [^, projection of low-mass excitations onto reso- 



nance states and diffractive inelastic interactions 39 



Corrections to the 2-string-excitation scheme arise at lower energies, if 
the excited masses which are randomly chosen according to eq. (|ip are below 
the cut-off values set for string fragmentation. If both outgoing states are 
below string excitation threshold, two hadrons (or resonances) are formed in 
the out-state 

h + h2^ hi + hi . (2) 

Production of two baryon resonances is the dominating process in pp colli- 
sions at the AGS energies of 10 to 15 GeV. Furthermore, even at the highest 
beam energies strings are actually rarely produced in nonannihilating events 
during the rescattering stage of nucleus-nucleus collisions. Therefore pro- 
duction (and absorption) of resonances in 2— i>2 processes effectively replaces 
the 2-string excitation component and fills up the inelastic cross section, in 
addition to s channel resonance formation and meson exchange processes in 
the t channel (see section |]). 

Let us assume now that above condition leading to a reaction of the 
type in (0) is met. The two produced excitations have to be projected onto 
hadronic states. Resonances are propagated explicitly in RQMD and may 
scatter themselves. All lightcone momenta which enter into eq. are cal- 
culated as if the ingoing hadrons would be groundstate hadrons, keeping the 
CMS energy the same. Thus the class of out-channels in (^ is populated 
with the same probability and cross section (given the same total cross sec- 
tion), irrespectively whether hi and /12 are excited states or not. The need to 
respect detailed balance requires some modifications which were neglected in 
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earlier versions of the RQMD model [Q. The improved model for such 2^2 
transitions will be described in the following. The flavor of each hadron is 
kept the same, an assumption which can be dropped eventually in further 
refinements of the model to allow for flavor-exchange between the collision 
partners. The out-channel is chosen randomly with proper weights to give 
the following cross sections for 2^2 where each term in the sum corre- 
sponds to the production cross section of a particular 2-hadron (resonance) 
state: 

a(fc,/^'2') = ao{s)J2pl (3) 



{2S^ + l)i2Sj + l)-SF{kl,ij) . 

i—l label hadron states, Pab denotes the CMS momentum and Sa the spin of 
hadron a. (To{s) is determined by the normalization condition. Depending on 
the initial choice - a nondiffractive or a diffractive inelastic interaction - the 
sum runs over all allowed states i and j or is restricted by the constraint that 
either hadron i or j is a groundstate hadron. SF{kl, ij) denotes a symmetry 
factor which is unequal to 1 only for k ^ I and i = j. In this case it is 1 /2 to 
respect the detailed balance relation. Since eq. @ is implemented in RQMD 
employing the Monte Carlo technique, the transition kl ii is assigned first 
a probability which ignores the symmetry factor. Afterwards the transition 
is accepted with probability 1/2. In case of rejection an elastic collision is 
realized. The physics content of eq. is rather simple. The out-states are 
chosen according to phasespace and available spin degrees of freedom, its 
functional form in accordance with the detailed balance relations. The basic 
idea is that statistics governs the population of out-states if many states are 
available in binary collisions. 

Extending the RQMD approach from hh to hA and AA collisions by al- 
lowing for multiple string excitations is straightforward. Each nondiffractive 
constituent quark interaction generates a new string if the energy is suffi- 



cient |jT3[. The cross section of interacting constituent quarks is given from 
the additive quark model ag^g = l/Q-cXiVAr which keeps the original projectile 
cross section constant while the projectile is traversing the target. Multiple 
interactions can be pictured as continuing 'undressing' of a constituent quark 
propagating through the target. 

It is an important question for the degree of baryon stopping and the 
achievable baryon densities in nucleus-nucleus reactions how the baryon num- 
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ber is shifted in multiple nucleon collisions with a target. The minimal- 
stopping approach is to concentrate all momentum which is not used up 
for target string excitation into a leading diquark which fragments like in a 
A^A^ collision. Historically, this was the first approach which has been ap- 



plied in string fragmentation models [^0], It is well-known by now that 
the minimal-stopping approach is at variance with the experimental data 



in ultrarelativistic pA and AA collisions |^. The shift of baryon number 
in rapidity is stronger. It contains a component which goes with the num- 
ber of additional projectile collisions. In the beginning several people have 
parametrized the dependence of the nuclear stopping power as a function of 



the number of additional collisions |^3|. A model for baryon stopping which 
contains a strong dependence on the number of collisions was suggested in 



33] and is applied in RQMD. It is in good agreement with available baryon 



measurements of the NA34, NA35 and NA44 group in 200AGeV collisions 
. In this approach the light cone momentum fraction of the baryon after 



fragmentation is determined from the number of spectator constituent quarks 
which it keeps in its wave function. The topology of string excitations in mul- 
tiple baryon collisions and the distribution of the ingoing momentum between 
the active degrees of freedom (strings and constituent quarks) are displayed 
graphically in fig. |[ Each further inelastic interaction in a target removes 
an additional constituent quark from the spectator remnant replacing it by 
a sea quark with zero momentum (in this approximation). It may even hap- 
pen that none of the original constituent quarks ends up in the outgoing 
baryon. In this case a partonic sea diquark moves at the forward end of 
one of the projectile strings, and the baryon momentum is determined from 
string fragmentation. 

The ingoing hadrons are complex objects themselves. Multiple collisions 
of a hadron at high energy are not sequential collisions of a single object, but 
in each collision a different hadron component is involved. Simple considera- 
tion of time scales tells that multiple soft interactions of the same object are 
highly suppressed. A fast particle which needs on the order of 1 fm/c in its 
rest system to finish the first interaction has left the target far behind which 
suppresses a subsequent interaction. This argument can be given a more 
rigorous meaning by analysing planar diagrams with multiple 'ladder' ex- 
change . Using the same line of argument multiple quasi-elastic collisions 
(2 2) and inelastic diffractive interactions should be suppressed as well at 
high energies. Of course, the basic reason behind the suppression of multiple 
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collisions within short time intervals is of quantum-mechanical nature. The 
uncertainty relation associates a finite time and distance interval to an in- 
teraction which can be characterized by exchange of energy and momentum. 
Recently a first step was done to take this effect into account in RQMD. 
Since the particle propagation in RQMD is realized semi-classically, only an 
approximate solution can be found in this framework. It was discussed a 
long time ago by Low and Gottfried that a classical space-time concept of 



propagation can be applied for fast particles because rapidity and lon- 
gitudinal position are commuting in this limit. Therefore one is allowed to 
specify t and y after a collision. The space-time interval after which a 
quasi-elastic or diffractive collision is considered as finished is now defined by 

^min = ^Coll + , (4) 

assuming that the particle is moving in forward (=+) direction. (A corre- 
sponding relation holds for the collision partner.) If the collision cannot be 
finished before the next collision is going to take place, the first collision 
being 'too soft' is discarded. The collision point xcou is defined from the 
minimum distance value which the two classical hadron trajectories can have 
in their 2-body center of mass frame. Aa;"*" is calculated using 

Aa;+ ~ l/Ap- , (5) 

with Ap~ the absolute value of the lightcone momentum which the hadron 
has picked up from the collision partner. (As usual in RQMD the longitu- 
dinal component is defined in the 2-body CMS by the direction of particle 
motion before collision.) Relation is motivated by the Heisenberg relation 
Ax^ ■ Ap^ > 1. The proportionality factor is presently set to 1. Applying 
the criterion of relation to all 2— i>2 collisions which are generated by 
RQMD in the dynamical evolution of nuclear collisions suppresses soft col- 
lisions. By their very definition, not many particles are produced in such 
collisions. Therefore the resulting effect on the final particle yiels is small 
in A(200AGeV) on A collisions, on the few percent level. The quantum- 
mechanical suppression of very soft multiple collisions in AA interactions 
at other beam energies is currently studied. The results will be presented 
elsewhere. 
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2.2 Color ropes in RQMD 

The quarks which have been polarized from the sea in course of the initial 
interactions are receding from each other. Their masses are usually small 
compared to their lightcone momenta and neglected. Therefore these quarks 
are moving on the lightcone and form the sources of the chromoelectric field 
filling the region between the forward and backward moving charges. The 
field is compressed into tubes with the tube cross section kept independent 
from the field strength. Thus a collection of tubes is spanned in AA collisions. 

First, I am going to describe the physics of a single chromoelectric flux- 
tube. The formation of tubes with stronger than the elementary string fields 
could be of relevance for AA collisions but also for soft multiparticle pro- 
duction in hadron-hadron collisions at collider energies |4^. The created 



SU(3)-valued color fields inside each tube are added coherently. The ele- 
mentary triplet charges of quarks and antiquarks in one of the receding rope 
'condensator' plates are coupled stochastically to the total SU3 color charge 
as the source of the rope field: 

(p,g)®(l,0)^(p+l,g)©(p-l,g + l)©(p,g-l) , (6) 

with the statistical weight given from the dimension of each multiplet 

d{p,q) = Y{p+l)-{q + l)-{p + q + 2) . 

The possible couplings of an arbitrary SU^ charge with an elementary charge 
are graphically displayed in fig. |^. 

The chromoelectric field of the rope is determined in the fiux-tube picture 
by Gauss' law: 

E"{p,q)-A{p,q)=gF{p,qr • (7) 

{p, q) characterizes the multiplet of the QCD charge. F(p, a = 1, 8 are the 
generators of the SU(3) group in the corresponding representation. It follows 
that the energy density of a rope field is proportional to the eigenvalue of 
the Casimir operator C{p,q) = F(p, g)"" ■ F(p, In RQMD the transverse 
size of a rope A{p, q) is being kept independent on the representation of the 
source. This is required in the fiux-tube model to get scaling of the rope 
tension with the Casimir eigenvalue which has been observed for both SU2 



and SU^ gauge groups g^l • Thus the rope tension, the energy per unitlength. 



also scales with the Casimir eigenvalue 

Kip,q)=3/A-C{p,q)-Kei . (8) 
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The tension of a triplet flux-tube (elementary string) Hei can be related to the 
Regge slope parameter a' which gives 0.9 GeV/fm. It is well-known that the 
flux-tube picture emerges from QCD with static quarks in lowest order of the 
strong coupling expansion Recently the assumed flux-tube poperties ~ 
dominance of the longitudinal electrical field components, scaling of the field 
strength with C {p, q) and charge independence of the transverse extension of 
the field - have been confirmed by a calculation of Wilson loops for higher 
dimensional charges in 3-dimensional SU2 lattice gauge theory |49]. Note 
that these new lattice results provide important information concerning the 
dynamics underlying confinement. If confinement is due to a bulk property 
of the QCD vacuum (like as the pressure in bag models), then the transverse 
size of higher dimensional flux-tubes is expected to increase. Consequently, 
the rope tension would scale less rapidly, as the square root of the Casimir 
eigenvalue only. Thus phenomenological models of confinement like the bag 
model are not compatible with the recent results of lattice simulations. 

Quark-antiquark pairs are created from the chromoelectric field and screen 
the original field [^. In a first approximation the total pair creation rate 
and the flavor composition can be calculated employing Schwinger's vaccuum 
persistence rate for a constant electric field . This is done by treating the 
pair creation as a tunneling process in the semi-classical WKB approximation 
(see Appendix). Employing a local density approximation, the calculated 
pair creation rate determines stochastically the space-time points in which 
the field strength is degraded by created quarks. This procedure general- 
izes the Artru-Mennessier scheme |^ to the situation of fields with variable 
strength. Of course, some modifications of the decay probability are expected 
for very short times and near to the sources of the rope field [^| . The pair 
creation probability per unit time and length is obtained by integrating over 
the transverse rope area vr-r^^^g with transverse rope radius set to 0.8 fm. The 
radius parameter is fixed by the requirement to get reasonable production 
rates from elementary string decays in comparison to e^e~ hadroproduction 
data. Strange quark production is easily enhanced by increasing the field 
strength, because the mass difference between the light flavors becomes irrel- 
evant for sufficiently large flelds. In contrast, 3-diquark creation inside the 
rope stays rather weak in the average. The small probabilities to create a 
diquark pair by the rope fleld have their root in the assumed 2-step process 
as outlined in Ref. [0]. 

The collective fleld is gradually degraded, because the initial 'macro- 
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scopic' spatial separation of charge and anticharge cannot be maintained. 
The first mechanism to degrade the field strength is quark pair production. 
Pair creation points, the created flavor and transverse momenta are sampled 
stochastically. The local field strength is calculated self-consistently, taking 
screening of the original field by already produced quark pairs into account. 
The rope tension is lowered in the forward lightcone of each break point 
from k{p, q) to g'). E.g., (p', q') = {p, q — l) if a quark is pulled out of the 
vacuum and attracted by the charge (p, q) . 

There are two other processes which degrade the original field strength, 
turning points of quarks in the rope end plate and crossing points of two 
color charges inside the rope. All these processes are displayed together in 
fig. ^. The turning point of a quark is determined from the condition that 
its original momentum has been used up by propagating under the force of 
the rope. Its momentum loss per unit-time is given from 

dp ^ ^ K{p,q) 

dt p + q ^ ' 

with the sign depending on the direction of motion. Of course, this is the 

direct generalization of the situation in string decay to account for the finite 

total momentum. Using eq. (|^) independent fragmentation is recovered in 

the limit that the fragmentation products of two strings have no overlap in 

rapidity space. | In this case the initial quark momenta show complete 

mismatch, e.g. one string is specified by {P~^,6~), the other by (5+,P~) with 

5"'"'~^P+'^. The region with non-triplet field strength shrinks to zero for 

5+'^^0 and P^~—>-oo. Similarly, those parts of string world sheets which do 

not overlap with other strings can make no contribution to build up a region 

of larger field strength. In RQMD these nonoverlapping parts of generated 

strings are split off from the beginning in the rope fragmentation process. 

The reason is the following. As displayed in fig. ^ two quarks (anti-quarks) 

may coalesce into a diquark in the 3 (3)-representation during the process of 

rope charging. Guided from the general coalescence picture one might expect 

that such a diquark would break if the momentum mismatch is too large 

between the two constituents, and an additional meson would be created. I 

have checked by introducing a 'reasonable' parameter for diquark break-up 

^ The hadron rapidity distribution from string fragmentation is approximately constant 
within some Umits ymm,max- The limits are related to the energy- momentum of the string 
via the relations ymin,max~^'0- P '^ -ln(2m), with m a typical hadron mass. 
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that practically the same results are achieved from rope fragmentation as in 
the parameter-free default procedure. 

The third process of field degradation, crossing trajectories of two color 
charges which do not form a singlet, has no analogue in an elementary string 
decay. It had also not been considered in early studies of the materialization 
time after which the rope field is degraded to zero [B3|. The two crossing 
charges need not form a white state as long as there are other color charges 
available which can neutralize their charge (see fig. Of course, two triplet 
or two antitriplet charges always cross without being able to form a color 
singlet. The probability that quark and antiquark form a white state is 
given in the model as {1+p' ■ with p' {q') being the number of 3 (3) 

charges available to neutralize the antiquark (quark) charge after crossing. 
This probability becomes very small for regions of high field strength. Thus 
color is confined globally in a rope, but not locally. In contrast, color charges 
are neutralized always locally in a fragmenting string. 

The quark pairs are produced in a rope with zero longitudinal momentum, 
but afterwards they are accelerated in the force field of the outer charges (see 
fig. The accelerating force is given by the difference between the rope ten- 
sions before and after pair production k-k', the same force which is driving 
the pair creation process. The classical trajectories are calculated neglect- 
ing the finite quark masses as it is usually done for string fragmentation. 
Therefore all charges are moving with the velocity of light all the time. This 
simplifies the calculation considerably and is, in fact, the only lorentzinvari- 
ant propagation without reference to the global rope rest system. 

Finally, all original and newly produced quarks will end up in a color 
singlet with a corresponding partner. The color singlets are projected onto 
the basic hadron multiplets, with the same relative weights as for string 
fragmentation. The generation of three-quark systems and projection on 
baryon states require some care Independent choice of three quark 

flavors in the rope endplates tends to overpopulate baryons with quarks of 
unequal flavor. Thus each chosen configuration is assigned a proper weight 
to avoid an unphysical fiavor SUq breaking. The positions of the breakpoints 
are slightly readjusted that all hadron momenta which are determined from 
the momentum sum of their constituents fulfill the mass shell constraints. 
The whole rope fragmentation scheme is constructed in a way that total net 
fiavor and energy-momentum are conserved. The hadron formation points 
are calculated from the quark trajectories. A formation point is defined 
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defined by 

x^=x^{Yo-Yo)-p^/{2k,i) , (10) 

where x^(Yo-Yo) is the first meeting point of the two quark constituents 
forming a hadron with momentum p^. In case of a triplet charges as the 
source of the rope field this prescription reduces to the corresponding defini- 
tion of formation points in string decay (cf. fig. |I]). 

It is clear from relations (|^) that in the process of rope charge formation 
a triplet and an antitriplet charge may eventually form a color singlet and 
do not contribute to the total rope charge (see also fig. 0). In the most 
extreme case the total charge is zero and the region of the rope is field-free! 
The probability that this may happen is rather small, however |Q. The 
quarks which form a color singlet in the rope endplates are projected onto 
hadron states of the basic fiavor SU^ multiplets. The momentum of such a 
hadron is given from the sum of the quark light cone momenta (if the hadrons 
are assumed massless, in the real calculation up to a small correction). The 
hadron formation point cannot be determined from the rope dynamics to 
which the quark constituents do not contribute in this approximation. It 
is usually expected that formation point and momentum are related in soft 
production processes, approximately ~ P^, which is used here with the 
scale factor k^i- 

The fragmentation of a single color rope as it is implemented in the 
RQMD model has just been described. In this picture several fiux-tubes 
are formed in very energetic AA reactions which may cover the whole trans- 
verse area of overlap between the two ingoing nuclei. How is the fusion of 
strings into several ropes realized in RQMD? At first the strings are gener- 
ated independently. The hadronization of each string is calculated in a Monte 
Carlo-type fashion, and the information is stored. This decay is considered as 
'virtual' until a first 'would-be' hadron of a decayed string is emerging at its 
formation point. All other strings which have been generated so far are now 
examined whether the transverse coordinates of their origin have a smaller 
distance to this first string than the fiux-tube radius rmbe- 'Transverse' and 
'longitudinal' refers to the coordinates in the rest system of this string. Thus 
the specification of the fusion process is lorentz invariant. Of course, the 

^ This process was neglected in the first RQMD calculations which included rope for- 
mation pT[ |. It changes the results from rope fragmentation on the order of 5 %, because 
the probability to form a color singlet by statistical combination of 2 randomly chosen 
(anti-)triplet charges is 1/18. 
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longitudinal direction coincides usually approximately with the beam direc- 
tion. Furthermore, it is checked whether the world sheets which are swept 
out by each QQ pair during half a period of string motion overlap in the 
t-z plane. Note that this criterion is trivially satisfied at asymptotic energy, 
because all string origins will be identical in this case. If no other string is 
found which fulfills the criteria the string decay is considered as 'real' and all 
hadrons emerging from its fragmentation are propagated subsequently. All 
strings which have been accepted are collected, and their independent decay 
is rejected. Instead, they fuse into a rope. 

In principle, the interactions between fiux-tubes can be more complicated 
than in the model adopted in RQMD, e.g. transverse 'communication' be- 
tween different tubes. So far such interactions are not taken into account. 
It should be noted, however, that the recent lattice results provide some 
justification for neglecting these interactions in a first approximation. Trot- 
tier and Woloshyn have shown that the color ropes in lattice gauge theory 
do not expand in transverse dimensions. This result has been obtained in 
the static limit, with an infinite amount of time available and with only the 
nonperturbative vacuum outside the tube. 

The need to go beyond independent string fragmentation for an under- 
standing of the early stage in ultrarelativistic AA collisions is generally ac- 
cepted nowadays. Of course, the appropriate replacement by more sophis- 
ticated approaches in line with QCD is currently under debate. It should 
be added that the string fusion approach of Refs. p2|, ^ based on the dual 
parton model is similar to the concept which is realized in RQMD. A major 
difference is that in the DPM calculations only fusion of at most two strings 
is considered so far. Fused strings break only as a whole in these DPM Monte 
Carlos while in RQMD stepwise field degradation is introduced. The RQMD 
approach seems better suited in case that many strings overlap and fuse. 
Nevertheless, possible effects of a different fragmentation procedure could be 
compensated partially modifying the fusion strength. A direct comparison 
of results from the two string fusion approaches which will be undertaken 
lateron shows rather similar trends for central Pb on Pb collisions at CERN- 
SPS energy. The beam energy may be too low that single-particle observables 
display sizable sensitivity to details of string fusion based models. 



16 



3 hh interactions in the rescattering stage 



The interactions in the hadronic resonance gas stage are described by binary 
coUisions between hadrons. Although preequihbrium processes are present in 



AA colhsions , the interactions are usuallly of nonasymptotic type, in the 
energy region of resonance production, formation and absorption. Most of 
the work in the RQMD framework has been devoted to develop a reasonable 
model for all so-called nonexotic reactions. In these reactions an s channel 
resonance can be formed as an intermediate state. A basic motivation is 
that the interaction strength is much larger here. (The energy-momentum 
entering the denominator of the corresponding T matrix element can be near 
to a pole.) Experimental justification comes from large differences of cross 
sections in nonexotic versus exotic channels, e.g. KN versus KN or tttc with 
isospin or 1 as compared to isospin 2, at low and medium energies. If 
ingoing hadrons form a state with quantum numbers not allowed for a sin- 
gle hadron by the quark model (exotic state), the low and medium energy 
interaction is solely given by t channel hadron (Reggeon) exchange. Only 
nonexotic reactions do allow quark-antiquark annihilation in the entrance 
channel. This process is very important for the flavor dynamics in AA colli- 
sions. If the system starts with strangeness below chemical equilibrium val- 
ues, quark-antiquark annihilation in nonexotic reactions drives the system 
towards chemical equilibration. Therefore this section presents a detailed 
description of the RQMD model for the nonexotic reactions. A modeling 
of these processes is clearly needed, because hadronic resonances which are 
treated as quasi-particles may interact themselves. For strangeness creation 
these processes are actually much more important than the interactions solely 
between groundstate hadrons. It was a very early observation in hadron 
physics that bumps show up in the energy dependence of cross sections which 
can be identified with the excitation of discrete Breit-Wigner type resonance 
states, e.g. the A(1232). Hadronic interactions at somewhat larger energies 
are dominated by quark exchange and annihilation (describable as Reggeon 
exchange) whose energy dependence are given by a negative power of s. At 
even higher energies Pomeron exchange in elastic interactions - exhibiting 



^ For instance, ttN interactions are insufficient to explain the A and K enhancement 
in S induced coUisions at 200AGeV, in contrast to e.g. the conjecture in Ref. The 
strangeness suppression in nN colhsions at the relevant invariant mass above f .6 GeV/c 
is rather similar to pp collisions at 200AGeV. 
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some kind of universality - induces an approximately energy independent 
total cross section for hadron-hadron collisions. These three components - 
with their characteristic and different energy dependences - can be identified 
with three pieces of interaction in RQMD: formation of discrete resonance 
states, annihilation into a resonance or string 'continuum' and excitation of 
two (or more) strings which was described in the section before. The im- 
portance of each component varies for different hadronic reactions. Four 
classes of interaction can be distinguished according to the number of ingo- 
ing (anti-)baryons {B denoting baryon, M denoting meson): BB, BM, MM 
and BB. The charge reversal invariance of strong interactions can be used 
to generate all other interactions of antibaryons. MM and MB interactions 
in states with nonexotic quantum numbers get some contribution from each 
component. In contrast, BB annihilation may occur only in the continuum, 
because the minimum invariant mass is too large. 

A description of the interactions which are induced by quark exchange 
and annihilation has to interpolate smoothly from the low energy region of s 
channel resonance formation to the high energy interactions which are better 
described in terms of t channel Reggeon exchange. Of course, the old dis- 
cussion about 'duality' of these two descriptions has never lead to satisfying 
results ^6[. For the construction of the hh interaction in RQMD related 
to quark annihilation and exchange I follow a pragmatic approach, adding 
together a few lowlying resonances and a Regge-type parametrization. The 
formation of s channel resonances is calculated in the MB and the MM 
sector from multichannel Breit-Wigner formulae. It is assumed for RQMD 
that resonance formation determines completely the interactions in nonex- 
otic channels up to some CMS energy ^/so■ This statement is not completely 
true as will be explained lateron. (There are some small corrections due to 
t channel background processes.) Above Sq when a description in terms of 
discrete resonance levels becomes invalid - the 'annihilation' cross section 
(for vrvr, vrA^, KN, etc.) is assumed to decrease with energy as Since 
the tails of the Breit-Wigner resonances decay much faster with energy, the 
resulting gap to the total annihilation cross section (o"^„„) is filled by anni- 
hilation into a 'continuum' which serves as the bridge to Reggeon exchange, 
at high energy supposedly dominant over s channel resonance formation. Q 

^ The cross section related to Reggeon exchange is dubbed here 'annihilation' cross 
section for terminological convenience, though part of the physics may really be quark 
'exchange'. I shall come back to this point later. 
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The choice of the exponent -1/2 for the power of s in the energy de- 
pendence of (TAnn IS motivated by the existence of energy thresholds in 
2 — i> 3,4,... reactions. They can be described by multiple Reggeon ex- 
change and effectively increase the maximum value -1 for the power of 
s in 2— i>2 Reggeon exchange processes as determined from the g trajectory 
(and, of course, from experimental data). One can test the assumed energy 
dependence by looking at reactions which are tied to Reggeon exchange and 
not to high energy processes determined by Pomeron cuts. For instance, 
the energy behaviour of the inclusive cross section K~p — > A + X - one of 
the rare cases in which flavor tagging allows for rather good distinction - is 
consistent with in the relevant energy region [s > 4GeV^). 

Since the annihilation cross section decreases with energy, a gap to the 
total cross section (atot) opens up. atot is either given from experimental 
data or calculated using the additive quark model [^. The physics filling 
the gap is the high energy component of hh interactions, Pomeron exchange 
for high energy elastic interactions and Pomeron cuts in inelastic collisions 



3J] . It follows from the considerations above that the energy dependence of 
the high energy component in reactions with measured total cross sections 



(e.g. TcN, KN) turns out to be approximately proportional to (l-yso/s). 
This functional dependence on energy is used to switch on the high energy 
component of the total and elastic AQM cross sections for reactions with 
unmeasured cross sections. 

After the general idea has been presented how the annihilation cross sec- 
tions are constructed as a sum of Breit-Wigner and continuum contributions, 
I am going to discuss the various classes of reactions specifically. The general 
scheme is realized somewhat differently in different channels depending on 
the level of experimental knowledge. 



3.1 Formation of s channel resonances 

Let us turn first to meson-meson interactions. Each of the three vector 
mesons p(770), _ft'*(892) and 0(1020) decays into two pseudoscalar mesons. 
It is experimentally well-known that p(770) and i^'*(892) dominate the phase- 
shifts in the p wave for vrvr, 1=1 and vri^, 1=1/2 scattering and can be de- 
scribed well by isolated Breit-Wigner resonances above a small background. 
The formation of these 3 and 26 additional meson resonances in the mass 
region up to 1800 MeV/c^ is taken into account in RQMD. Note that not 
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all SU3 flavor nonets are complete due to lack of experimental information 
about resonance masses and widths. However, all lowlying multiplets - the 
scalar, the two axial vector and the tensor mesons - are completely included. 
The groundstate mesons and the resonances of these multiplets are propa- 
gated explicitly in RQMD, while the other resonances do appear only in the 
intermediate states. 



Godfrey's and Isgur's quark model calculation |]5^ which is rather suc- 
cessful in explaining measured branching ratios of resonance decays has been 
taken to extract resonance couplings to unmeasured decay channels. Partial 
decay widths of resonances which can decay only off-shell into a 0, have been 
related to measured branching ratios employing flavor SU^ symmetry and a 
correction for phasespace kinematics. This applies to the K* resonances of 
the lowlying tensor and two axial vector nonets. 

In RQMD Breit-Wigner type multi-channel cross sections for s channel 
resonance formation are usually summed up incoherently. Exceptions are tttt, 
1=0 and 7tK, 1=1/2 reactions in the s wave. Here interferences are taken 
into account, because very broad resonances are present in these channels 
(/o(1400) and i^rQ(1430)) which interfere with a strong background from at- 
tractive interactions [0]. In addition, the narrow /o(975) state interferes 
destructively in the scalar-isoscalar channel with the other /q resonance and 
a background. The phases of resonances and background are added in the 
elastic channels 

So = s, + J2^R , (11) 

R 

which is the appropriate way to get an unitary S matrix as long as only 
one channel is open. The background phase 6b can be calculated from one- 
meson exchange in the t channel employing the K matrix formalism (see 
the subsection |37^ ). = arctan is the phase attributed to each 

resonance. In particular, unitarity is important for the vrvr interaction in the 
scalar-isoscalar channel. Just below invariant mass of 1 GeV the phaseshift 
goes through 180 degrees due to the strong destructive interference effects. If 
more than one channel opens up with increasing energy (in vrvr interactions 
the KK channel), the situation becomes even more involved. In order to 
ensure continuity of the S matrix across the thresholds eq. ([TT|) is also used 
above particle production threshold. The inelasticity parameter 77 in the 
diagonal elements of the 5* matrix S^^ = rji exp (2i(5j) is calculated from the 



20 



unitarity condition 

T.si:si = i (12) 

k 

with the (small) non-diagonal elements given from addition of Breit-Wigner- 
type T-matrix elements. No relative phase between resonances and back- 
ground is introduced which would allow the construction of a completely 
unitary multi-channel 5* matrix. This is possible in the framework of the 



Davies-Baranger formalism . The results in the scheme adopted here and 
in this more complicated unitarization approach are rather similar, however 

In all meson-baryon interactions with quantum numbers allowed for a 
baryonic state from flavor SU3 resonances may be formed. So far, however, 
Q* formation is neglected in RQMD. The properties of nonstrange baryon 
resonances - with isospin 1/2 (A^*) and 3/2 (A*) - have been experimentally 
well explored in ttN collisions, and of hyperon resonances (with isospin 
and 1) in KN interactions. The knowledge about multiply strange baryon 
resonances (S* and Q*) is poor, however. A generalized Breit-Wigner formula 
for the cross section is used to calculate the s channel resonance formation 
probabilities (for A^*, A* and S*): 

^.(. ^ ^ (2Ji;+l) TnjMB) ■ r^(tot) 

■ ■ V (25i + 1){2S, + 1) ■ (yi - mn)^ + r«(tot)V4 " 

(13) 

b{s) is a normalization factor which renormalizes - for A^* and A* for- 
mation - the Breit-Wigner sum to a given absolute cross section. (For S* 
formation incoherent addition gives b{s)=l.) The renormalization is done 
for nonstrange baryon resonances, because the total cross sections in vrA^ re- 
actions are known. The total resonance formation cross section calculated 
with eq. ([I3[) in vrA^ is actually the total vrA^ cross section minus an elastic 
background (up to 5 mb) determined by a consistency condition that the 
sum of background and elastic resonance decay equals the measured elastic 
ttN cross section. TjilMB) denotes a partial decay width into the channel 
with meson M and baryon B. The decay widths are -y/s-dependent via the 
relative momentum in the CMS of the decaying resonance p 

TniMB) ~ niR/^s ■ p(2'+i)/(i + o.2 ■ {p/p{mji)f) . (14) 

This ensures a correct threshold and high energy behaviour of a particular 
channel. The particular resonance which is formed in a meson baryon anni- 
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hilation is chosen randomly with a weight given by the corresponding term 
in the sum of eq. ([T3|). 

In the A^* and the A* channel the sum in eq. (|T3|) runs over all resonances 
with a mass below 2 GeV/c^. The resonance masses and decay parameters 
as implemented in RQMD are in accordance with the listing of the Particle 
Data Group It has been checked that an application to vrA^ collisions 
reproduces approximately the hyperon-kaon production cross sections and 
the total pion yields in inelastic collisions ||3^ . 

In the S* channel formation cross sections for 6 discrete resonance states 
are incoherently added with couplings to meson-baryon states given from 
S'f/3 flavor symmetry. (Note that flavor SU^ is broken by different M and B 
masses in each octet.) The F and D parameters specifying the strength of 
symmetric and antisymmetric coupling are taken from the literature . 



Resonance formation as given in eq. ([T3|) is not realized for hyperon reso- 
nance formation (A* and S*). Instead, the experimentally measured exclusive 
cross sections for charge exchange {K~p K^n), elastic interactions and hy- 
peron production {KN <-> ttA/ttE) have been tabellized in RQMD. The S* 
resonance parameters are not determined very well. Furthermore, at low en- 
ergy the KN interaction is rather complicated, e.g. due to the existence of 
interfering resonances. 

The lifetime of formed resonances is stochastically chosen according to 
an exponential decay law. The average lifetime (without collisions) is given 
from the inverse decay width. The angular distributions in decays of formed 
resonances are certainly an area in which more work is needed. All transport 
models for AA collisions so far had ignored this problem, which is already 
present in the 'prototype' reaction t^N —>■ A (1232) -kN , by decaying the 
resonances isotropically. While the problem can be easily solved for isolated 
resonances, it becomes more involved for overlapping resonances. The dif- 
ferent partial waves should be summed coherently which would generate a 
diffraction-type peak at small scattering angles. As a first step isotropic res- 
onance decay is supplemented in RQMD by a second component, a Gaussian 
Pt distribution. This ensures that the average pt which is generated does 
not exceed the standard value of 400 MeV which is usually assigned to a 
produced quark pair in soft production processes at high energy. 



^ Some of the resonance parameters have been changed shghtly as compared to the 
values given in because they were refitted together with the modified continuum 
component to -kN data. 
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3.2 Transition from s channel resonances to Reggeon 
exchange 

The formation of discrete s channel resonances, which represents the low 
energy component of quark- ant iquark annihilation in hh interactions, has 
just been described in some detail. Such a description is expected to break 
down in the intermediate CMS energy region s = 4 — 16 GeV^. Resonance 
level densities are getting too high and resonance widths too large that a 
formulation in terms of discrete resonance excitations could be meaningful. 
Furthermore, due to these difficulties the empirical information about dis- 
crete states gets poorer with increasing invariant resonance mass. On the 
other side, Regge theory is spectacularly successful in this energy region. I 
do not discuss here the speculation that an infinite tower of s channel reso- 
nances could effectively generate the same interaction as t channel Reggeon 
exchange (see e.g. the Veneziano model []63|). 

The main emphasis is put in this subsection on a discussion of 5=1 chan- 
nels. The continuum component in MB reactions is of utmost importance 



for nucleus-nucleus collisions due to preequilibrium processes |T^, while it is 
suppressed kinematically in MM collisions. For instance, here the processes 
are above threshold for ss creation. There is one more reason to model 
carefully MB interactions in this energy region. Strange (anti-)baryons are 
a very promising signature of collectivity in ultrarelativistic nucleus-nucleus 
reactions. After they are produced in an AA collision, they have to propagate 
eventually through a cloud of nonstrange mesons fl^ . The strength of tran- 



sition rates like Ap KN or EK KY is determined by the probability 
of moving the s quark from the baryon to the meson side (and vice versa). 

The total strength and the energy dependence of the Reggeon exchange 
component follow from the ansatz which was introduced above: 

<7Regge{s) = (cTtotis) - aBw{s) " ^(s)) for S < Sq (15) 

= (o-toi(s) - CTBwis) - At(s)) ■ \Jso/s S > Sq 

for ttN, KN and vrS interactions. Eq . ([15|) is applied for MB states of good 
fiavor quantum numbers (isospin, hypercharge). Interference effects of am- 
plitudes in different states are neglected. The ^/so parameter is set to 1.8 
GeV in the N*, 2 GeV in the A* and Y*, and 2.15 GeV in the E* channel. 
At(s) is given from the elastic background cross sections for vrA^ or the sum 
of the tabehzed cross sections KN — > KN, vrA, vrS for KN collisions. Ac- 
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cording to eq. (pTSf ) the contribution of the Breit-Wigner resonances in the 
N* / A* /E channels have to be subtracted to get the total 'Reggeon exchange 
component'. Note that the tails of the A^* and A* resonances above sq are 
calculated from eq. (|^) with the b{s) values 'frozen' to the value of b at s=so. 

No experimental data are available for the total Svr cross section. Its value 
above sq is taken from the additive quark model (AQM). Lack of experimental 
information is the reason that fewer S* resonances are included in the Breit- 
Wigner sum of eq. (plsf ) than in other channels. This is 'corrected' in eq. (|T5|) 



by multiplying the constant AQM value by yso/s for s smaller than sq down 
to the S*(1530) resonance, if the Breit-Wigner sum gives a smaller value for 
the cross section than this parametrization. 

There is no Reggeon exchange component in the A^* and A* channel for 
energies below s = sq. (The annihilation cross section is filled up by s channel 



resonance formation according to eq. (|T3D.) In the A* and S* channels the 



Reggeon exchange component is present down to lower energies, because the 
difference of the total cross section to the sum of the tabelized cross sections 
becomes nonzero, even below invariant mass of 2 GeV/c^. 

The total Regge exchange component as given by eq. (^) can be subdi- 
vided into single Reggeon exchange (2 — > 2) and into multi-Reggeon ex- 
change (2 — y 3,...) diagrams. Multi-Reggeon exchange diagrams with 
n > 2-body final state may become important to account for some frac- 
tion of multiparticle production at higher energies. However, such processes 
are suppressed in the hadronic rescattering stage of AA collisions Their 
presence would be indicated by an additional component to the produced 
particle yields nonlinearly increasing with the number of participant nucle- 
ons for which there is also no hint from experimental data. Before the relative 
weights for these two classes of processes in RQMD are discussed, it is use- 
ful to describe how a particular 2-body final state in a reaction of the type 
MiBi — > MoBo is selected out of the sample of all possible states. 

The exclusive cross sections of the Regge exchange component 2 — ^2 have 
been given the following functional form: 

.iM.B. ^ MM - ^.(,,^^^,^)(,,^^^,) E^K^)- (16) 
||Cmb,/|' (l - exp(-(pMB/^)')) {{2Sm + + l))r(i?)j 
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Eq. ([16|) contains several factors - in their ordering from the left to the 
right - related to incoming flux and spin averaging, absolute normalization 
(a/(s)), isospin projection (Clebsch-Gordon-coefficients squared), threshold 
phase space factor, internal spin degrees of freedom, suppression parameter 
r{R) of strangeness creation or quark exchange between M and B, and a 
suppression for exchange of Reggeon R as determined by the intercept of the 
corresponding Regge trajectory aR. If the absolute normalization is ignored 
for the moment, the relative probability of some final state is determined 
rather simply. Up to some statistical factors which express the 'matching' of 
the isospins and the available spin degrees of freedom there are only three 
factors which contain dynamics: the phase space factor, the strangeness pa- 
rameter r{R), and the Regge factor. Since the Regge factor tends to 1 at low 
relative CMS momenta p and p', the approach assumes kind of 'universal' (or 
average) coupling of a meson-baryon 2-body state to the intermediate system 
with baryonic quantum numbers. The phase space factor cuts down on the 
production probability near threshold, but then also goes quickly to 1. Its 
functional form is motivated by the available phase space for Schwinger type 
Pt creation in flux-tube decays which is Gaussian distributed, with width 
(T=400MeV. The parameter r{R) is assigned a value different from 1 only to 
account for ss suppression (0.3 in a connected and 0.09 in a disconnected 
quark-line diagram) and to suppress exchange of quarks between M and B 
in some processes (0.5 for exchange of baryon number or strangeness). The 
first parameter is the usual strangeness suppression which is characteristic 
for soft particle production. The last parameter value reflects the immobil- 
ity of the heavier strange quark and of the leading spectator diquark. The 
model presented here is very similar in its spirit to Vandermeulen's statisti- 
cal approach |Q, and Mundigl, Vicente Vacas and Weise's doorway model 
[|65|] , both for meson production in low energy pp annihilation. This is not 
really a surprise, because the common denominator is the overwhelmingly 
large number of open states which precludes any kind of detailed dynamical 
treatment. 

Some additional strangeness suppression is 'hidden' in the Regge factor 
which is smaller for strange Reggeon exchange. It is a general observation 
that the strangeness suppression in nonexotic hh — >hh interactions is some- 
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what stronger - in addition to phase space suppression - than 'asymptoti- 
cally' in soft hadronic physics {P{s)/P{u) is ~ 0.3 for multiparticle produc- 
tion). See e.g. the suppression parameters for kaon channels in pp annihi- 



lation [^4|, ^ which are needed to fit the experimental data. A stronger 
ss suppression is expected from microscopic calculations which give a larger 
suppression for heavy mass quarks if a color flux-tube has a small length 



or breaks very early |5^. However, it cannot be ruled out that part of the 



stronger suppression effect comes from the neglect of quark exchange dia- 
grams. Their inclusion - in reactions without strange quark in the initial 
state - favors nonstrange final states. It has been found for MB interactions 
that an average strangeness suppression parameter of 0.2 can reasonably well 



describe total strangeness creation in vrA^ reactions |]53|. Therefore not too 
much freedom is left for a pure quark exchange component in such reactions. 

The Reggeon exchange suppression factor contains the momentum pji as 
a scale at which Reggeon exchange becomes important. Thus the parameter 
choice reflects the interpolation from the s channel resonance region to the 
high energy Regge approach. Aaji is defined as 

AttR := 2 - 2aR - ao . (17) 

For each MiBi — >MoBo reaction only a single Regge trajectory is taken into 
account, the one with largest intercept. The vertices which are used in the 
RQMD calculations are presented in tables |1] and |^, together with the pa- 
rameters for the Reggeon intercepts. 

The multi-Reggeon component is represented in the model by a hh an- 
nihilation into a color string which decays into more than two hadrons (res- 
onances) according to the standard string fragmentation scheme. In prin- 
ciple, the probabilility whether a 2-body final state (2^2 process) or an 
n-body final state {n > 2) is chosen in a Reggeon exchange process is deter- 
mined by the normalization factors aj{s) in eq. ([I6|) . However, aj{s) is not 
specified explicitly in RQMD. Instead, the correct Regge behaviour for 2^2 
processes at large energies provides an implicit normalization of the total 
2^2 Reggeon-exchange cross section. Regge theory requires all 2^2 cross 
sections to behave (approximately) as 

at large energies. With the parameter choice ao=l the cross sections of 
eq. (p6|) show this behaviour under the assumption that a/(s) approaches a 
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constant value at high energies. 

The correct asymptotic behaviour for 2—>-2 transitions is reahzed in the 
following way. For invariant masses of the colliding MiBi pair below some 
value y/si=2.8 GeV (plus 0.15 GeV for each s quark in the B* state due 
to flavor symmetry breaking) it is assumed that all final states are 2-body. 
Above Si the 2-body final state is suppressed by an additional power 1/2 

(^ReggeC^ ^ 2) < aRegge{s) ■ \J Si/ S , (18) 

and the appearing gap to the total Reggeon-exchange cross section is filled up 
by choosing an iV>2-body final state. The < sign in eq. ([l^) is to be under- 
stood in the sense that crijegge(2— *-2) is assigned at first a value which is equal 
to the right hand side of eq. (0). The normalization according to eqs. (p!5|) 
and (0) would give asymptotically a 1/s dependence for aReggei/^^'^) if 
the equal sign in eq. (^) would hold. This is actually the correct energy 
dependence if Reggeons of the q trajectory with the largest intercept are ex- 
changed (the Ug ~ 0.5). On the other side, cross sections for other transitions 
MiBi — > MoBo may drop faster with energy due to a smaller intercept pa- 
rameter aR of the exchanged Reggeon. In order to avoid that cross sections 
characterized by a comparably larger intercept aR simply blow up due to 
the normalization condition in eq. (|TB|), a fraction of the total 2^2 Reggeon- 
exchange cross section goes actually into many-body final states. This will 
be explained in the following. 

After the decision to choose a 2-body state has been made in RQMD, a 
particular state has to be selected, respecting its relative weight according 
to eq. (|16D. This is realized in a Monte Carlo- type fashion using the rejec- 
tion method. After some randomly chosen out-state has passed all other 
acceptance tests, the Reggeon factor 

PRegge{s) = {l+Pp/p%) ^""^ (19) 

is taken as a probability for final acceptance of the 2-body state. Below 
Si, Monte Carlo type rejection of a 2-body channel is followed by choosing 
another 2-body state until at last an out-state passes all acceptance tests 
successfully. Above Si, the acceptance probability is again given by PRegge{s). 
However, the probability to examine another 2-body state is frozen at the 
value l-pRegge{,si). As a third possibility, an A^>2-body final state may be 
chosen, with probability 

probability(iV > 2) = pRegge{si) - PRegge{s) for S> Si . 
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This prescription guarantees the continuity of all cross sections across si and 
the correct energy dependence of all 2—^2 cross sections asymptotically. 

Which hadrons are included in the in- and out-states of eq. (|16|)? All 
nonstrange baryon resonances with a pole mass below 2 GeV/c^ whose ex- 
istence is established by the Particle Data Group can be produced. All 
baryon multiplets with these A^* and A* states have been completed by 
including strange baryon resonances either according to the data tables, 
flavor SUs relations or by applying a simple formula for resonance mass 
(m = m{N*/A*) + ■ 150MeV) and width (Tr = dRT{N*/A*) with (Ir 
equal 0.85 for S*, 0.62 for S* and 0.45 for Q*). Thus an unwanted flavor 
SU^ breaking due to the poorer experimental knowledge of resonances in the 
strange resonance sector is avoided. In a second step the formed resonance 
states decay subsequently. The branching ratios of those resonances whose 
properties are not taken from the data tables or flavor SU^ relations have 
to be specified. They are determined similarly as given for the intermediate 
doorway state in eq. (0), however, without Reggeon factor: 

T{B* MB) ~ (20) 
\Cmb,i\' (l - exp(-{pMB/cry)) {{2Sm + 1)(2^b + l))r(i?*) . 

Here the parameter r{B*) is related to strangeness again. Its value is set 
to 0.2 if ss is created in the decay (1 otherwise). Of course, strangeness 
creation in a second step is much more improbable than in the first stage of 
interaction given by eq. (pIBI) . 

With the combined model of s channel resonance formation and t chan- 
nel Reggeon exchange one can describe reasonably well on the order of 60 
measured exclusive meson-baryon reactions (vrA^ and KN), including their 
energy dependence. Usually the results agree with exclusive 2 and 3-body 
data better than within a factor of two. Inclusive yields, especially of strange 
hadrons (kaons, hyperons, etc.), are described even much better (within 10 
percent). The constructed model does not aim for a 'perfect' fit to experi- 
mental data. Instead, the intention was to introduce only a minimum number 
of parameters whose values have some physical meaning like strangeness sup- 
pression, the scale which separates the resonance from the Regge region and 
so on. Therefore the model can be meaningfully extrapolated to ingoing hh 
combinations with experimentally unexplored properties of their interactions. 
No new parameters have to introduced if the n or K and the nucleon are re- 
placed by other meson-baryon combinations, in particular resonances. Their 
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low and medium-energy interaction are completely specified by the same 
model, essentially given by eqs. (|I3|) and (p^Gf). The interactions fulfill the 
property of detailed balance. Detailed balance follows from the time reversal 
invariance of strong interactions. A cross section aio for a MiBi — > MqBo 
transition is related to the cross section in the reverse direction by 

p\2SMi + msBi + i) 

p'\2SMo + mSso + lf''' ■ ^ ^ 

Thus the model avoids the pitfalls of equating cross sections with ingoing 
resonances to groundstate hadrons as used in early RQMD calculations 
and in other transport models 

After the discussion about the quark annihilation processes in the con- 
tinuum has focused so far on B=l channels, I would like to discuss shortly 
its role in channels with total baryon number zero. Such reactions can be 
either BE annihilations or MM interactions. The MM annihilation into 
a continuum is constructed very similarly to the approach outlined above 
for MB interactions. The continuum annihilation starts to open up at 
so=((l-6+ns-0.15) GeV)^, because the Breit-Wigner resonance tails decay 
faster than s^^^^ which again is the assumed energy dependence of the to- 
tal annihilation cross section. Note that the coupling between 2-body in- 
and out-states is not reggeized yet. Instead, the couplings are constants 
with strangeness annihilation and creation suppressed by the factor 0.2 as in 



3^, A parametrization of the pp annihilation cross section is used which 
was fitted to experimental data ^6|. The construction of the other BB ab- 
sorption cross sections is described in [Q. In principle, one should put the 
same effort into a detailed description of the out-states as discussed for MB 
interactions. However, the choice of the outgoing states in the continuum re- 
gion is of minor importance for the reaction dynamics of heavy ion collisions. 
MM collisions around an invariant mass of 2 GeV/c^ are rather rare in A A 
collisions |13|, because too much internal kinetic energy is required. Clearly 
BB collisions are not kinematically suppressed in this invariant mass region. 
Therefore the final B yield is very sensitive to the annihilation strength. 
However, it matters only whether an antibaryon survives in baryonic matter 
or not. The specific choice for the particle production model in BB annihi- 
lations has not much infiuence on the meson dynamics itself. Note that the 
typical B/M ratio in central AA reactions at 200 AGeV is on the one percent 
level only. So far the RQMD procedure for the choice of the outgoing states 
in the B = system is to form a mesonic string which decays according to 
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the standard rules of string fragmentation. In addition, it is enforced that G 
parity is conserved in nonstrange meson transitions by rejecting out-states 
with wrong G parity. 

Particle production at energies in the 'continuum' region, i.e. beyond the 
regime of identified resonance decays, can be naturally viewed as a tunneling 
process of quark pairs leading to subsequent breaking of a color flux-tube 
Thus the annihilation processes in the continuum are translated in RQMD 
into the language of (Yo-Yo type) string decay. This determines the space- 
time structure, most importantly the formation points of outgoing hadronic 
states. 



3.3 Meson exchange processes in the t channel 

Here I am going to discuss the component in RQMD which describes t chan- 
nel exchange driven processes in MM interactions with two mesons or res- 
onances in the final state. Meson-exchange processes in MM interactions 
are generated by a flavor SU3 symmetric Lagrangian which couples a vector 
meson to two pseudoscalar mesons. Q 

The empirical knowledge about meson-meson interactions is rather poor 
due to the inherent difficulty to prepare a mesonic 'target'. Only vrvr and 
TcK interactions are well explored. It has been outlined above that formation 
of the p and K* resonance are dominating the interactions in these chan- 
nels. The decay of the vector meson resonances into two pseudoscalars is 
approximately described by a flavor SU^ symmetric interaction term in the 
Lagrangian 

Lint = -{i/2)GvTT{[P,^^P]V^) = 2Gvf^JkP^^^PJVJ!; . (22) 

P and V denote the pseudoscalar and vector meson matrix [|7I|]. Flavor 
symmetry is broken only by the mass differences between mesons within the 
multiplets. 



In addition to vector meson decay, the SU^ invariant Lagrangian in eq. |22 
generates transition rates for reactions PP PP, PV PV and PP * 



^ In addition, RQMD includes parametrizations of cross sections involving nucleons 
and the A(1232) resonance which were calculated using the One-pion-exchange model 
(OPE) e.g. NN iVA, or OPE-inspired parametrizations of some measured 

flavor creating processes like NN A{1232)YK 0. These processes are of not 
much importance in 200AGeV collisions, the topic of this paper. 
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VV by exchanging a meson in the t channel. Gy is fitted to the p decay width 
and set to 3.02. All relativistic Born diagrams for fiavor changing processes 
and for all important exotic channels (which have no other interactions at low 
energy) have been calculated and implemented in RQMD. The interaction 
vertices are supplemented with monopole form factors 

which express the finite size of the interacting mesons. The wrong energy de- 
pendence in reactions with vector meson exchange is 'cured' by multiplying 
each PPV- vertex with (s/sq)^"""^^/^ for s > Sq {y/s^=lGeV+mi+m2). SU3 
invariant values for the cutoff-parameter A are chosen, A= 3.0 GeV^ (1.5 
GeV"^) for attraction (repulsion) in the s wave. The calculated interactions 
give remarkably good agreement with the measured phase shifts in tttt and 
irK interactions ||S^,[^-[|73|. As an example the s wave phaseshift is shown 



for TTTT scattering with 1=0 and 1=2 (cf. fig. Note that the 1=2 channel is 
exotic, while the 1=0 channel contains additional contributions from s chan- 
nel resonances, in particular the /o(975) and the /o(1400). A completely 
unitarizing scheme - eq. (0) below KK threshold and its generalization, 
the Davies-Baranger formalism, above threshold - has been applied to calcu- 
late the combined effect of t channel background and s channel resonances. 
The real-valued matrix elements for one-meson-exchange have been identi- 
fied with the K matrix which has been decomposed into their components 
from different partial waves. The background S matrix for partial wave / is 
constructed from the K matrix via: 

Si^ = {l + iK'){l-iK^r' . (24) 

While a standard parametrization is used for the /o(1400) decay width, the 
KK molecule picture is employed for the /o(975) [^. Furthermore, the 
/o(975) decay width into KK is analytically continued below threshold like 
in Ref. |^ in order to generate the sharp cusp observable in the scalar- 
isoscalar phase shift just below y/s=l GeV. 
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4 Flavor production in central Pb(160AGeV) 
on Pb collisions 

The two basic building blocks of interactions in RQMD - color ropes and 
hadronic rescattering - have been presented in the preceeding sections. In 
this section the consequences of these collective interactions for the flavor 
dynamics in the most central collisions of two lead nuclei at a beam energy 
of 160AGeV will be discussed. The impact parameters for the calculated 
Pb(160AGeV) on Pb collisions have been selected between and 1 fm, rather 
central collisions. The calculations with the RQMD computer code have been 
done in different modes, the default mode and two alternative modes. The 
first alternative mode is defined by switching off all collective interactions, 
no rescattering and no rope formation ('NN mode'). The second non-default 
mode is without rescattering but rope formation included. 

4.1 Particle multiplicities 

The results of the different RQMD calculations for the final hadron yields 
are presented in table |^. All members of the basic pseudoscalar meson nonet 
and (anti-)baryon octet have been kept stable here with the exception of 
the rl which decays already during the dynamical evolution generated by 
RQMD. Some extreme scenarios are a scaling of the totally produced hadron 
yield with either the number of participants or the number of binary colli- 
sions. These scenarios lead to an dependence of the produced particle 
multiplicities with a values of 1 and 4/3. One would get the a- value 4/3 
from RQMD in the 'NN mode' by neglecting finite-energy effects. In this 
case the final particle yield would simply scale with the number of binary 
collisions, because in each collisions new strings can be created and each 
string decay gives asymptotically a constant rapidity density. The number 
of created quark pairs in the 'NN mode' of RQMD is 3061 which gives more 
than 7.2 created hadrons per participating nucleon. This can be compared 
to the corresponding numbers in elementary "p-p collisions which is 4.9 (at the 
slightly higher beam energy of 200AGeV). In terms of an parametrization 
this means an a- value of 1.07 by comparing Pb+Pb and reactions. This 
value is still far away from the 'upper limit' 4/3. Obviously the finite-energy 
effects - finite string masses and mutual deceleration of projectile and target 
- change the naive estimate considerably. 
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How do the production rates of quark pairs change if rope formation is 
taken into account in Pb(160AGeV) on Pb colhsions? The formation of co- 
herent chromoelectric fields which fill the space between the receeding color 
charges sets different initial conditions for particle creation than the incoher- 
ent superposition of elementary strings. The effect of ropes on the conversion 
rate of field energy into qq pairs, the resulting rope lifetime and particle mul- 
tiplicities have been discussed already in the literature [H, Neglecting 
partially or completely the screening effect of already produced charges and 
employing various approximations, in particular boost invariance, simple re- 
lations for these variables as a function of the rope field strength have been 
derived. The authors of Ref. @] derive a relation valid in their model 

^ 2.12 K^oo , (25) 

with F the meson multiplicity and its argument the so-called foldness. The 
foldness K is defined by the number of (anti-)triplet charges which screen a 
rope field completely. K=p+q for a rope whose source is a color SU^ charge 
in the representation. Eq. ( pS] ) implies a strong suppression of particle 
production ~ compared to production from N independent strings. 

Note that random charging without any constraint leads to the average K 
value 

The multiplicities from rope fragmentation in RQMD cannot be calculated 
analytically, but show similarly a trend that particle production is strongly 
dampened with increasing rope charge. The assumption that each hadroniz- 
ing string breaks just once seems more appropriate for an estimate of the 
rope effect at 160AGeV on multiplicities than the opposite extreme of an 
infinite number of break points which is a consequence of assumed boost 
invariance. Under this assumption the number of created quark pairs is de- 
pleted if ropes are formed, because the minimum number of quark pairs to 
screen the original field is smaller (a/ZV versus A^). Thus the total number of 
quark pairs including the quarks of the source charge is given simply as 2A^ 
from string decays, +A^ from rope fragmentation. 

Neglecting finite-energy effects grossly overestimates the influence which 
rope formation can have on the whole dynamical evolution of the system. 
After the hadrons which contain ingoing constituent quarks have been sub- 
tracted from the average number of primary hadrons in elementary hh col- 
lisions at IGOGeV is around 2 per collision partner. Subsequent collisions in 
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a nuclear target produce effectively even fewer hadrons per collision which 
come from smaller strings or from excited resonance decays. 

The surface and the finite thickness of the ingoing nuclei play a nonnegli- 
gible role as well for rope formation. The Lorentz contracted lead diameter in 
the CM. frame is approximately 1.4 fm. The corresponding passage time for 
the collision partner has therefore a value which is larger than the hadroniza- 
tion scale (1 fm/c). Strings or ropes which were formed early in the collision 
may have already hadronized before the projectile has completed its passage 
through the target. All these effects tend to suppress the importance of rope 
formation and favor the other sources of secondaries, elementary strings and 
resonances. 

The RQMD approach to rope formation and hadronization is superior 
to schematic calculations, also because it respects the simple constraints 
arising from nuclear geometry and finite beam energy. The yield of produced 
secondaries decreases from 3006 to 2571 in the calculated central Pb on Pb 
collisions by allowing string fusion into ropes. This 15%-effect on the total 
multiplicity from rope formation is tiny as compared to the suppression effect 
implied by eq. It will be interesting to study the energy dependence 

of the multiplicities in case of rope formation. The dampening of particle 
production due to formation of strong chromoelectric fields should become 
much more pronounced at collider energies. 

Inclusion of hadronic rescattering leaves the total produced particle mul- 
tiplicities practically unchanged (2512 in comparison to 2571 without rescat- 
tering). These values correspond to an a- value of 1.04 which is even closer to 
1 yet than in the 'NN mode' of RQMD. The reason is that particle number 
conserving processes - in RQMD 2^ 2 - dominate the rescattering stage. 
Particle conservation in the expansion stage of AA collisions has been conjec- 
tured by other authors earlier based on general grounds (kinetic equilibration 
together with G parity conservation which forbids to change the pion number 
by one unit in collisions of nonstrange mesons) |7^ and is confirmed by the 
microscopic transport calculations. The small decrease of the particle mul- 
tiplicity due to rescattering is connected with the strangeness enhancement 
(see below) and the related smaller feed-down of the lightest hadronic degree 
of freedom, the pions. One should keep in mind that the strange hadrons. A, 
Ks, etc. are considered stable here. 
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4.2 Rapidity distributions 

Here I shall focus on the shapes of the rapidity distributions for different 
particle species. The rapidity distributions for negatively charged hadrons 
calculated in the different RQMD modes are displayed in fig. 0. It is clearly 
visible that the change in produced particle multiplicities is concentrated in 
the central region. The collective effects leave practically no tracks in the 
fragmentation region of projectile and target (|?/cms|>2). At these rapidities 
the three rapidity distributions - from the 'NN mode', with ropes and with 
additional rescattering fall on top of each other. This can be understood 
from the discussion in section p^.2| . The rope fields influence mainly the par- 
ticle multiplicities at rapidities which are covered by many strings. Also the 
small absorption effect which rescattering has on the particle multiplicities 
is concentrated around midrapidity. On the left hand side of fig. the net 
proton (p-p) rapidity distribution is shown for comparison. The p-p rapidity 
distributions are remarkably similar in the three RQMD modes. Rescatter- 
ing makes the final net baryon distribution only slightly narrower. The small 
effect of hadronic rescattering on the baryon stopping in the present RQMD 
calculations differs markedly from results in the string model QGSM (with- 
out string fusion but rescattering included) |7^. The effect in the calculation 
presented here is also somewhat smaller than in early RQMD studies [BC . 
The difference is mainly caused by the more realistic angular distribution 
in baryon resonance decays (see section ^Tl] for a discussion). Isotropic de- 



cay as employed in Ref. tends to shift baryons rather rapidly into the 
central region. One can see from fig. |^ that the p-p rapidity distribution is 
significantly broader than the produced particle distribution. However, the 
baryon stopping mechanism in RQMD is much more pronounced than in 



other approaches [^, |79[. In RQMD it shows no minimum at midrapidity 
but a plateau in the central region. 

The comparison of the rapidity distributions of net protons and negatives 
in fig. 1^ demonstrates that the single fireball picture is inconsistent with the 



calculated rapidity distributions, like in the smaller S+A reactions |T^. The 
primordial nucleon rapidity density would have a width between 0.44 and 
0.52 units at freeze-out temperature between 140 and 200 MeV. The heavier 
a particle, the narrower is its rapidity distribution in a thermal fireball. The 
reverse ordering found from RQMD is indicative for the presence of longitu- 
dinal flow in the system. One could imagine different mechanisms driving the 
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longitudinal flow, internal pressure (as in Landau scenario) or a transparency 
effect (as in Bjorken scenario). RQMD in which the secondaries are initially 
created by the boost-invariant decay of a longitudinally stretched tube is con- 
structed very similarly to a Bjorken-type approach. Some differences arise 
from flnite-energy effects which are clearly visible from flg. |^. There is no 
plateau in the produced particle rapidity distribution at all. Furthermore, 
the Bjorken ansatz of identifying space-time and momentum-space rapidity 
is 'softened' in string models resulting in a flnite width of the local rapidity 
distribution [Q. The calculated rapidity distribution of negatively charged 
hadrons is narrower in central Pb on Pb than in S(200AGeV) on S colli- 
sions. The width of the distribution is 1.6 units of rapidity versus 1.9 which 
makes a difference of 8% if they are normalized to the initial rapidity gap 
between projectile and target. The narrowing arises in the model from the 
stronger attenuation of the ingoing nuclear constituents in the heavy system. 
It is interesting to note that Landau-type hydrodynamics with all matter 
initially at rest would show the opposite effect. Fig. || displays the calculated 
rapidity distributions of strange baryons and mesons. The distributions of 
strange baryons are markedly narrower than the proton density, 1.7 versus 
2.1. This cannot be attributed to the mass differences which results only in 
small width differences in a thermal picture (< 0.05). Taking longitudinal 
flow into account does not explain this effect either. In RQMD it reflects 
differences in the production dynamics. The probability that a baryon car- 
ries flnally a strange quark is correlated with its rapidity. Baryons at central 
rapidity have suffered simply more collisions, initially with nucleons from the 
other nucleus, lateron with secondaries. Spatial inhomogeneity of the flavor 
composition is no real surprise, even in collision of truely heavy ions. On the 
other side, parametrizations of 'hydrodynamical flow' under the homogeneity 



assumption have been quite popular for studies of light ion reactions |81 



One can take the RQMD calculations presented here as an indication that 
results in such simple models have to be taken with some grain of salt. 

The flnal antikaon and kaon distributions have a slightly smaller width 
(1.3-1-4) than the hyperons as a result of the generated longitudinal flow. In 
contrast, the antibaryon rapidity densities have much smaller widths (1.1-1.2) 
than mesons and baryons as can be learnt from the results presented in flg. ^ 
Again, this effect is incompatible with a homogeneous thermal flreball - with 
or without longitudinal flow - and reflects different production mechanisms 
in the model. Since antibaryons are mostly produced in rope fragmentation 
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(see below), they are primarily formed in position and rapidity space where 
many strings overlap. Baryon pairs emerge only from the regions of highest 
energy density. In contrast, due to smaller production threshold ss pairs 
are created also in the diluter regions as well. It had been already noted by 
the experimentalists in the presentation of first data for S(200AGeV) on S 
collisions that the different shapes of A and A dN/dy distributions provide 
evidence for differences in the production mechanisms |^3[. The systematics 
of measured strangeness and antibaryon production in S on A collisions at 
beam energy of 200AGeV [51| can be exploited to understand better the 
dependence on initial baryon and energy density. The latter variables are 
related to the final baryon and meson rapidity densities. It was argued in 
[ P^ that the narrow antibaryon distributions in S+A reactions f\ arise as a 
convolution of two effects, string or equivalently energy density, which is the 
determining factor in the model for the achievable rope field strength, and 
the smoothly increasing baryon density (in position and momentum space) 
towards the target rapidity region. In contrast, in central Pb on Pb reactions 
the antibaryon rapidity densities get slightly broadened due to absorption, 
because the rather flat net baryon density for Pb on Pb is centered around 
midrapidity. 



4.3 Strangeness enhancement 

How large is the fraction of created strange quarks relatively to the light fla- 
vors from RQMD? The fraction of created ss pairs normalized to the average 
of uu and dd pairs 

2N{ss) 
" N{uu) + N{dd) 

turns out as 11.3% in the 'NN mode' which should be compared to 10% in pp 
collisions at 200GeV, a small increase. f\ The rate of produced kaons remains 
practically unchanged in the 'NN mode' if they are normalized to the pion 
yield. For instance, the /tt ratio is 6 percent (with isospin averaged pion 
yield) while it is 5.6 % in pp collisions. Several competing effects may give 
rise to small deviations from the pp result. Because of energy degradation the 
average string masses are lower in multiple collisions of a projectile hadron 

^ Note that the p distributions for the calculated reactions with O and S projectiles in 
fig. 5 of Ref. jl^ are scaled by a factor of 10 (not 20 as indicated in the figure). 

^ It is be neglected here that rj and ?/ meson carry some amount of hidden strangeness. 
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impinging on a nuclear target. This effect is unfavorable for strange particle 
production. On the other side, the ingoing flavor has to be used up if only 
one string is excited in course of a nucleon collision. If several strings are 
attached to an ingoing interacting constituent quark, all of the additional 
strings have color charges from the sea, possibly an s-quark, at their end 
points. Furthermore, multi-step processes of the type B—* 5*— > B** which 
are most relevant near particle production threshold may still play a small 
role even at a beam energy of 160 AGeV. 

The strange baryon rate normalized to the total baryon number increases 
considerably already in the 'NN mode'. The ratio of A/B including all feed- 
down contributions is 0.05 in pp collisions from RQMD which is in accordance 



with experimental data |8^. It increases by nearly a factor of 2 to a value 
of 0.09 in central Pb on Pb reactions. This effect can be understood rather 
simply from the multiple collision effect on baryons (see the diagrams in 
fig. 0). While only 1 or 2 of the original valence quarks are replaced in 
elementary NN collisions, the large majority of valence quarks are being 
replaced in the outgoing baryon states by quarks from the sea in Pb on Pb 
collisions at small impact parameters. Thus the probability increases with 
the number of collisions which a baryon undergoes that it will carry a single 
or even multiple units of strangeness after the interaction. A similar effect 
is observed in the VENUS approach to AA collisions which includes the 
mechanism of 'double-string' formation . 



A strong strangeness enrichment compared to the 'NN mode' and the 
calculated values for pp collisions is expected by including the collective in- 
teractions into the calculations. The suppression of heavy quark production 
in strong chromoelectric fields is weakened. The different barrier penetra- 
tion factors for virtual quark tunneling are approaching the same values. 
Rescattering processes tend to enrich strangeness in nuclear collisions, be- 
cause - starting under conditions in which strangeness is undersaturated - 
nonstrange quark pairs are preferentially annihilated in resonance and string 
formation processes and sometimes replaced by a strange quark pair in the 
decay. Associated production plays the most important role for strangeness 
enhancement in the hadronic resonance gas. The difference between a A 
baryon and nucleon mass is just 178 MeV, but twice this value between 
anti-kaon and pion mass. Furthermore, the preequilibrium contribution to 
the secondary interactions is stronger in the meson-baryon than the meson- 
meson sector . This is qualitatively easily understood, because the baryon 
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rapidity distribution is broader than the meson distribution. This holds also 
locally which is responsible for the hard tail in the collision spectrum of 
baryons. Transport calculations with the RQMD model show good agree- 
ment with available data for reactions with light ion projectiles which are 
sensitive to this production mechanism (net-A rates, / K~ enhancement) 

mm- 

Most of the strangeness is carried in the final state by (anti-)kaons which 
makes their yield normalized to the total produced particle multiplicity to a 
signature for the strength of collective effects in AA collisions. It has been 
discussed already some time ago based on studies of S induced collisions at 
200AGeV that rope formation leaves the absolute yield of produced kaons 



practically unmodified |TT[. The situation does not change if collisions with 
heavy projectiles are considered (see table 0). The relative enhancement 
due to the weaker strangeness suppression is completely counterbalanced by 
the depletion of total multiplicity. The results in the DPM based string 
fusion approach SFMC confirm these conclusions and are surprisingly close 
to RQMD, even in terms of absolute numbers. The 47r yields for central Pb 
on Pb collisions given in Ref. are (with RQMD values in parantheses) : 
negatively charged hadrons without fusion 924.4 (964.3), but with fusion 
included only 806.2 (829.2). In contrast, the kaon yields change much less, 
K+ 90.1 versus 88.4 (79.2/79.0) and R- 71.5 versus 66.8 (53.0/50.4). 

The Rs value increases in RQMD from 11.3 % to 15 % if the strings 
can fuse into color ropes and even to 24% if the hadrons interact with each 
other after formation. One can learn from these numbers that most of the 
strangeness enrichment is produced in the hadronic stage and not during 
the quark matter evolution in central Pb on Pb collisions. Furthermore, the 
strangeness enrichment in comparison to pp interactions is even stronger in 
central collisions with Pb than with S projectiles. Both aspects are actually 
connected as will be discussed in the following. The strangeness enhance- 
ment factor for quark pairs which are created in the rope field is practically 
projectile-independent in comparing central S and Pb collisions, namely 50 
%. This is again a consequence of the finite energy and length scales which 
tend to suppress string excitations and longitudinal overlap of strings after 
the first 2-3 collisions of a colliding nucleon. Therefore the prehadronic pro- 
cesses are not responsible for the differences in the final states of reactions 
with S and Pb projectiles. 

There are two reasons why the interactions in the hadronic stage are more 
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effective in enricliing strangeness in reactions witli larger projectile masses. 
Tlie first reason is tlie increased stopping power wliicfi shifts relatively more 
nucleons of tlie target towards tlie central rapidity region. This allows them 
to participate more easily in interactions with secondaries which are concen- 
trated around midrapidity. The prominent role of baryons in increasing the 
strangeness content of the system can be directly read off from table |^. The 
final baryons contain 124 s valence quarks, which are 45 % of the s quarks 
created in total. A nuclear hypermatter state is created in these collisions 
according to the RQMD calculations. With respect to strangeness produc- 
tion, the dynamics of baryons is far from being a perturbation to the whole 
dynamics in AA reactions at the energy of 160AGeV. |^ 

The second reason is related to the time intervals which the system spends 
in these two 'phases'. |^ After the quark matter and hadronic preequilibrium 
stages which last about 3 (4) fm/c in S (Pb) induced reactions have been 
completed the system stays in approximate kinetic equilibrium for the rest of 
the interaction time. The temperature drops only slowly which is a property 
of the resonance gas with its many degrees of freedom. The hadrons and 
resonances interact in this stage until freeze-out, for approximately 11 fm/c 
in central Pb+Pb, but only for 3 fm/c in S+A reactions. Thus the weight of 
interactions in this third phase is more pronounced in Pb induced reactions 
than in light ion (e.g. O or S) collisions, even with heavy targets. In the 
latter case transverse expansion becomes rather effective after local kinetic 
equilibrium is reached, because the transverse area of the system is smaller. 

4.4 Baryon pair production 

A strong strange antibaryon enhancement has been experimentally observed 
for central S collisions on S and heavier targets [^, 0, Binary inter- 
actions in the rescattering stage tend to conserve particle numbers. This 
excludes them as an important contributor to the creation of additional an- 
tibaryons if the total interaction times are on the order of 4-6 fm/c only. 

^'^ Consideration of baryon degree of freedom is sometimes neglected in theoretical stud- 
ies of AA collisions which is based on the small B/M ratio in the final state at 160- 
200AGeV. The important role of baryons for strangeness production which is found from 
RQMD has repercussions for the production of other particles whose production is inhib- 
ited by thresholds, for instance direct photons at large transverse momenta. 

A detailed account of the time evolution of the created matter in central 
Pb(160AGeV)-|-Pb collisions will be given elsewhere 
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Such lifetimes of the interacting system in S+A colhsions are suggested from 
interferometry data ^] consistent with the RQMD calculations. The 
failure of Boltzmann-type hadron resonance gas dynamics to account for a 
sizable enhancement of strange antibaryons in central AA collisions points 
towards the importance of production mechanisms in denser stages of the 
reaction, notably in the prehadronic stage. This is supported by a simple 
argument. The energy needed to create a baryon pair is on the order of 2 
to 3 GeV. Soft production processes which are characterized by the scale 
Aqqd ^ Ifm therefore require energy densities of a few GeV/fm'^ in order to 
overcome the BB suppression in comparison to meson production. 

Color ropes are a strong source of additional baryon pairs in ultrarela- 
tivistic AA collisions |Tl|. There are two distinct mechanisms how a rope 
may generate diquark pairs in the triplet representation (3 — QQ). These are 
quark coalescence in the rope end plates and quark pair creation with color 
mismatch in screening the field. The two processes are graphically repre- 
sented in fig. 0. In the process of forming the rope charge elementary triplet 
charges of quarks and antiquarks in a receding rope 'condensator' plate cou- 
ple stochastically to the total SU3 color charge of the rope. This provides a 
strong source of 3 — QQ's and 3-QQ's by quark coalescence. The statistical 
weight that two quarks form a 3-QQ compared to a 6 — QQ is just 1:2. The 
(anti-)diquarks in triplet representation will combine with a corresponding 
anticharge to form an (anti-)baryon. 

It turns out that in nuclear collisions at 160-200AGeV the probability 
of quark coalescence in the charge at the end of the rope is typically much 
larger (>90%) than due to diquark creation inside the field. This is very 
much based on the energy gain by optimal screening in course of charge 
creation ~(C(p,q)-C(p,q-l)) and/or ~(C(p,q)-C(p-l,q)) as compared to the 
nonoptimal case ~(C(p,q)-C(p-l,q+l)+2/3). Note that this result reflects 
also the rather short average length of the fiux-tubes which are created at 
these energies. Giant fiux-tubes with infinite length would be completely 
dominated by creation processes inside. 

Indeed, baryon pair production is strongly increased in central Pb on 
Pb collisions if overlapping strings are fusing into ropes. The total num- 
ber of created baryon pairs in the 'NN mode' turns out to be 27.7. Ropes 
increase this number to 84.3. The SFMC calculations show qualitatively 
similar trends, an increase (of 2p+A) from 16.7 to 83.9 due to string fusion 



23]. However, the flavor dependence is different in the two approaches. In 
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RQMD the increase is strongest for strange antibaryons which is reflected 
for instance in the increase of the A/p ratio. The trend is the opposite in 
the SFMC results. The RQMD result arises from the 'constructive interfer- 
ence' of weaker strangeness suppression and additional diquark production 
mechanism if ropes are formed. 

Subsequent antibaryon absorption in baryonic matter brings the final an- 
tibaryon yield from RQMD down again to 21.8 in the final state. This means 
that the strong antibaryon enhancement (a factor of 3) from the preequilib- 
rium quark matter stage is more than 'eaten up' by the subsequent absorption 
in baryon-rich matter. Only a quarter of the initially produced antibaryons 
survives the interactions in the hadronic gas environment until freeze-out. 
Like in the case of strangeness enhancement a comparison of the yields in 
the different RQMD modes of operation demonstrates that the hadronic stage 
has much more impact on the flavor composition in central Pb-|-Pb than in 
S induced reactions [IT3. The strength of the antibaryon absorption shows 



some flavor dependence. While the number of p is even a factor of 2 below 
the results in the 'NN mode', the final yields of strange antibaryons are still 
considerably enhanced. For instance, the S yield increases by a factor of 7.7 
compared to the result in the 'NN mode', even a factor of 13.3 if antibaryons 
would not be absorbed in the hadronic stage. Since both antibaryon and 
strangeness production are suppressed in elementary hadronic interactions, 
the rope formation effect is particularly strong for hadrons which carry these 
flavors in combination. Furthermore, the flavor dependence of antibaryon 
absorption reflects that ss pairs are less frequent in the system and their 
annihilation probability is reduced, which mirrors the suppressed production 
probability (cf. section p.2|) . 

Is there some way to disentangle the antibaryon enhancement effect due 
to color rope formation or some other collective effect like QGP formation 



in the early stage of the reaction from strong antibaryon absorption in 
the later hadronic stage? At least, it can be expected that the rather extreme 
model for antibaryon interactions adopted presently in RQMD which assumes 
a free pp annihilation cross section can be tested rather well in heavy ion 
collisions. Since the annihilation cross section is strongly energy-dependent, 
preferentially low-momentum antibaryons are annihilated. This effect is even 



more pronounced in heavy ion collisions at lower beam energies |91]. The 



use of the free pp annihilation cross section is certainly debatable in 
particular at low relative momenta. Since this cross section corresponds 
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at low energy to particle interaction distances clearly larger than 1 fm and 
thus exceeds the inter-particle distance in the dense stage, medium effects are 
expected to be very important. The real part of the antibaryon self-energy in 
baryon-rich matter may get substantial values as well, because vector meson 
exchange leads to attraction in addition to attraction already from scalar 
exchange. It will have some influence on the shape of the particle spectra, 
although the magnitude of this effect is unclear pTI . 



5 Summary and Conclusions 

String fusion into color ropes and hadronic rescattering - collective interac- 
tions in the preequilibrium quark matter and hadronic resonance gas stage 
of ultrarelativistic nucleus-nucleus collisions - have been modeled in RQMD. 
The system created in central Pb(160AGeV) on Pb reactions is character- 
ized by strong longitudinal flow. The nucleon momentum distribution shows 
the strongest elongation along beam direction, but with a maximum still at 
midrapidity. The broadness arises from partial transparency (corona effect). 
The final antibaryon source is clearly more concentrated around midrapidity 
than their antiparticles, kaons and pions. Antibaryons are produced only 
in the region of highest energy density, while mesons and even strangeness 
are created also in the diluter regions as well. If these RQMD predictions 
are experimentally confirmed, strangeness and antibaryon enhancement in 
ultrarelativistic heavy ion collisions cannot be described by only one source 
with homogeneous flavor composition. The strangeness suppression factor 
defined by the ratio of created quark pairs ss/ {uu+dd) is strongly enhanced 
by a factor of 2.4 in comparison to pp results. Color rope formation increases 
the initially produced yield of antibaryons to 3 times the value in the 'NN 
mode', and even stronger if they carry strangeness. Only approximately one 
quarter of the produced antibaryons survives because of subsequent strong 
absorption in baryon-rich matter. The differences in the final particle com- 
position for Pb on Pb collisions to S induced reactions are attributed to the 
hadronic resonance gas stage which is baryon-richer and lasts longer in the 
heavy system. 
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Appendix: Quark pair creation in a rope 

The path integral for the tunneling of a virtual charge pair is given in the 



WKB approximation as |92 



P{Pt) = |e-'l^^'lf , (26) 

where Sd = Jo^'" Piiz) dz is the action for the path across the classically 
forbidden region. Therefore longitudinal momentum and action are purely 
imaginary in this region, initially: 

pf+p? + ^2^0 . (27) 

The force acting on a created charge is given from the amount by which the 
field energy per unit-length is lowered due to screening of the original source 



0: 



F, = {k-k') , (28) 

with K the rope tension before and the tension after pair creation. After 
each of the created charges has moved over a distance z the energy balance 
reads: 



2F,-z = 2^pi{z)+pl + m^ , (29) 

for a charge with constant mass \Sci\ = 7rm^/4(/t— k'). However, it is assumed 
here that the mass varies linearly with distance which is motivated by the 
expectation that quark masses are 'current' on short distance scales and 
'constituent' masses with respect to the nonperturbative confining force. 

m{z) = mo + P ■ z, (30) 
z < Am/p , 

with mo the current quark mass (m^ = m^ = 10 MeV, m^ = 160 MeV), 
Am=350 MeV the difference to the constituent mass, and (3=0.355 GeV/fm 
the 'speed' of quark dressing. The parameters are fitted to give 0.1 for sup- 
pression of diquarks in comparison to quark production and 0.29 for strange 
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quark suppression in elementary flux-tube decays which are the values fa- 
vored by experimental data. 

The absolute value of the action Sd for the tunneling of a quark with 
linear interpolation between current and constituent mass is given by 



\Sr. 



h 



(31) 



h = ITT 



2A 



M Fj 



Hi) 
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where the following abbreviations have been used: 



A = JF^ - p\ 



B 



mot ■ (3 



M = mot 



A 



(AM mot ^ 



Z2 



mot + Am 



Casher, Neuberger and Nussinov (CNN) have rederived the exact Schwinger 



result for the vacuum persistence probability in the WKB approximation 



with the additional benefit of getting out a transverse momentum distribu- 
tion for the created charges. Following their derivation the pair production 
probability per unit-time and unit- volume in a uniform Abelian field is given 

by 



7 

47^3 



CO „ 



flavors ^^=1 



n 



(32) 



7 denotes the degeneracy factor from the color degrees of freedom. In total 
there are 2 x 3=6 different (anti-)quark color states. However, only three 
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of them can lower the field strength by screening a given rope charge (see 
fig- p. 

A charge vector in the SU3 multiplet {p,q) - with p>q for definiteness - 
could combine with a created (anti-)quark charge to a SU3 charge state in 
the multiplets {p-l,q) or (p-l,g+l) to lower the field strength. One 

can see from fig. [TT] that the sum of original SU3 charge and a color-'up' 
quark do not have only a nonzero component in the but also in the 

(p-l,g+l) multiplet. However, the latter possibility is discarded, because it 
is energetically unfavorable. Therefore the factor 7 in the pair creation rate 
eq. (|32D representing the color degree of freedom is set to 1 in the three cases 
given above {q^O assumed). Screening of a color charge with g = can result 
in one state with charge {p-l,q) or in two states with charge (7=2 
in this case). The two configurations with screened charges in the {p,q-l), 
(p-l,g) multiplets will give 3-3 color singlet states together with the residual 
rope charge (optimal screening). The third configuration (p-l,g+l) consists 
of a 3-diquark on one side and an anti-diquark with opposite color charge on 
the other side (color mismatch) and possibly some other charge. 

The color of a diquark created by screening {p,q)—*{p-l,q+l) will get neu- 
tralized in a next step by an additionally produced quark pair. This is the 
generalization of the 2-step process for baryon pair production in an elemen- 
tary fiux-tube 1^ to the case of stronger chromoelectric fields. The model of 
diquarks adopted here resembles very much the picture which one gets from 
strong coupling QCD ^3||. A baryon state is represented in configuration 



space as a system of three quarks, each of them connected to a junction, 
which couples them to a color singlet with the help of the Levi-Civita tensor 
e (Y-shaped string). An effective diquark carrying 3-color charge would con- 
sist of two quarks and the junction which should be treated as an additional 
dynamical degree of freedom in principle. In a diquark-creating tunneling 
process a piece of elementary color flux is replaced by the e-j unction and two 
'legs', each having half of the elementary string tension (from energy conser- 
vation). Thus the force creating a diquark pair in an elementary flux-tube is 
weaker by a factor 1/2 than the one which produces an optimally screening 
quark. This explains naturally the dynamical suppression of baryon pairs as 
compared to meson production in elementary flux-tube decays. 

The produced quark pair acquires some transverse momentum in the tun- 
neling process whose distribution can be calculated from eq. (^). However, 
one should not expect that a tube is straight and parallel to the distance 
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vector between source and sink of the electric field on a small distance scale 
since it is subject to 'roughening'. Zero-point oscillations of a tube's normal 
modes may provide another source of transverse momentum. Kokoski and 
Isgur discuss the effect of roughening for the breaking of a flux-tube in the 



strong-coupling limit |^5|]. On the other side, the effect should become irrel- 
evant in the classical limit of infinite electrical field strength. Here I follow 
the pragmatic approach to add two uncorrelated components for a produced 
particle's transverse momentum, one from tunneling and another one from 
unresolved transverse excitations which give for the absolute value: 

Vtr — Vo^ Vtunn • (33) 

The po component is Gaussian distributed and fixed by the requirement that 
the total created transverse momentum in elementary fiux-tube breaking is 
(approximately) 400 MeV. 
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Table Captions: 



Table 1: 

Reggeon exchange in meson-baryon collisions: planar diagrams. The Reggeon 
which is exchanged in the various meson-baryon interactions describable by 
a planar quark-line diagram can be read off from the table. Values for inter- 
cept ur, strangeness and quark exchange suppression parameter r{R) and 
momentum parameter pji are given here which enter into eq. (0). m denotes 
here a meson with strangeness S=0, k with S=l, irrespectively which mul- 
tiplet it belongs to. (In contrast, a K denotes a kaon and K* a i^'*(892).) 
The multiplet of which an exchanged meson is a member is indicated by a 
V (vector) or PS (pseudoscalar meson nonet). Alternatively, the parame- 
ters are specified as some 'average' of the V and PS exchange parameters. 
Note that in the uppermost diagram the flow of isospin quantum numbers as 
specified and G-parity determine whether a or a PS meson is exchanged. 

Table 2: 

Reggeon exchange in meson-baryon collisions: exchange diagrams. Contrary 
to the Reggeon exchange diagrams in Table 1 the reactions tabelized here 
include a quark exchange between ingoing meson and baryon. The notation 
is the same as in Table 1. 

Table 3: 

The produced hadronic state in Pb(160AGeV) on Pb collisions with impact 
parameters b<l fm, calculated in three different operation modes of RQMD 
2.1: so-called NN mode with ropes and rescattering switched off in the right 
column, rope fragmenation included in the middle column and ropes and 
hadronic rescattering both included in left column (default mode). All mem- 
bers of the pseudoscalar meson nonet and the baryon octet have been kept 
stable, except the r]'{958) which is decaying already during the dynamical 
evoulution generated by RQMD. 
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Figure Captions: 



Figure 1: 

Decaying hadronic string in position space (scaled by Kei)- Four hadrons A-D 
are formed, their formation points indicated by dots. Thick hnes give the 
trajectories along which the hadrons, respectively in the beginning the orig- 
inal constituent quarks, are propagated, pj^, PyoYo^ ^"^^ ^^e forward 
light cone momenta for the leading (spectator) quark, for string excitation 
and for the interacting (backward) quark, p^^ is the backward light cone 
momentum which has been transferred from the target. 

Figure 2: 

Schematic diagrams for string excitations in multiple baryon collisions: The 
primary string excitation is shown at the top (a). Either the interacting 
quark collides two or more times (c) or the spectator diquark of interaction 
(a) interacts (b). Iteration of diquark interaction is depicted in diagram (d). 
Here all original valence quarks are completely stripped off from the baryon 
which is emerging from the fragmentation process. Constituent quarks are 
symbolized by full dots (interacting quarks on left, spectators on right side), 
sea quark pairs by open circles. The ingoing light cone momenta - forward 
Pi (/ for interaction) and ps {S for spectator) and backward p^ (provided by 
the target) - are distributed onto outgoing constituents and string excitation 
iPf-Pb) as indicated just below each diagram. 

Figure 3: 

Possible SU(3) multiplets which can be built by a combination of {p,q) states 
and elementary triplet (antitriplet) states, p is the number of columns with 
one line, q the number with two lines - (p=2, q=2 in this example). Thus p 
is the number of 'quark-like' charges, q of 'antiquark-like' charges. There is 
no restriction in forming the total rope charge. It is energetically allowed to 
create a charge in rope decay only if it screens the original rope charge, i.e. 
the resulting charge is lowered (right above, middle below and -p>q assumed 
- middle above). Thus only (at most) three out of six quark or antiquark 
states in color space can be created in the rope field. 

Figure 4: 
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Schematic picture of quark pair creation in a rope: quark trajectories are 
displayed in the t-2;-plane. The field strength is characterized by the {p,q) 
values of the source (which is acting here from the right side) . It is indicated 
in the figure how charge creation and crossing of quark trajectories decrease 
the field strength inside the rope. Quark and antiquark may form a color 
singlet which splits from the rope (symbolized here by two horizontal lines). 
Two examples of hadron formation which are displayed here demonstrate 
that color is not necessarily locally confined in a rope field. Neighboring 
quark Qd and antiquark Qc form a color singlet. Such a topology would be 
always enforced in elementary strings ((p,g) = (l,0)) by energy conservation. 
In contrast, Qa and Qc travel quite some distance before they combine into 
a color singlet. 

Figure 5: 

Three mechanisms of chromoelectric field degradation: quark pair creation 
(a), turning point of a quark in the end plate source (b), and crossing point 
of two quark trajectories (c). The field strength is characterized by a pair of 
numbers like {p, q) which characterizes the charge acting from the right side. 
(The corresponding anticharge which is the sink of the fiux on the left side 
is a member of the {q,p) multiplet.) Square brackets like characterize 
the charge moving on a particular trajectory. While p and q can get assigned 
arbitrarily large integer values, a charge which is denoted by [i,j] in the 
diagrams belongs either to a triplet ([1,0]) or an antitriplet ([0, 1]). 

Figure 6: 

Phase shifts 6q and 6^ in tttt scattering: Comparison between calculation and 
experimental data. The data are taken from Refs. [^, |7^ . 

Figure 7: 

Rapidity distributions of negatively charged hadrons and net protons in 
Pb(160AGeV) on Pb collisions with impact parameters b<l fm, calculated 
in three different operation modes of RQMD 2.1: ropes and rescattering 
switched off (dashed fine), rope fragmenation included (dotted line) and ropes 
and hadronic rescattering both included which is the default mode (straight 
line). The rapidity is calculated in the equal-speed-system of projectile and 
target. The negatively charged hadrons include feed-down from weak decays 
(except from Ks and (anti-)A), the net protons from all weakly and strongly 
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unstable baryons. 



Figure 8: 

Rapidity distributions of strange baryons and kaons in central Pb(160AGeV) 
on Pb collisions. The different histograms are related to the different RQMD 
operation modes as explained in the caption to fig. 7. Note that the A 
distribution does not contain any feed-down from H decay or S^-decay. The 
S-rapidity distribution is calculated by averaging over all isospin states. The 
amount of splitting in the yields of the three states is small, however (see 
Table 3). The S baryon distribution in the figure contains the sum of both 
charge states. Again, both states are populated approximately with equal 
strength. 

Figure 9: 

Rapidity distributions of anti-baryons (p. A, H) in central Pb(160AGeV) on 
Pb collisions. The different histograms are related to the different RQMD 
modes of operation as explained in the caption to fig. 7. The S distribu- 
tion contains the sum of both charge states. No feed-down from weak and 
electromagnetic decays is included. 

Figure 10: 

The two different production mechanisms of baryon pairs from rope fragmen- 
tation are depicted here schematically. An antidiquark which is part of the 
original source charge of the color flux field may combine with an antiquark 
created by the field itself (left side). It is assumed here that the original 
charges making up the total rope charge move along the light cone. Hav- 
ing lost its momentum by pulling out the chromoelectric field an end plate 
charge turns its direction and gets accelerated again. After combining with 
a corresponding anticharge into a colorless state it will split from the rope 
without further interaction. This process is visualized here in the t-z plane 
(with z the direction of the electric field). A flux tube may create as well a 
diquark-antidiquark pair in the color triplet configuration (right side). The 
scheme employed here is usually called two-step (or sometimes 'popcorn') 
production mechanism. It was first suggested in Ref. fSlI] for the case of 



baryon production from elementary flux tubes. The strength of the flux is 
indicated in the figure by a pair of numbers like (p,g) which defines the rep- 
resentation of the charge source. Note that the notation {QQ)z to highlight 
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the quark content of the color charge means a (1,0) (or triplet) charge. 
Figure 11: 

The three configurations in color hypercharge and isospin space for screening 
of a charge in the multiplet {p,q) - here the state |(p,q),t=(p+q)/2,t^=- 
t,i/=(p-q)/3) - by a created (anti-)quark charge. With p > q which is as- 
sumed here there is one screening antiquark color (a cross enclosed by a 
circle) and two quark colors (each represented by a cross). 
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Diagram 
(planar) 



Subclasses / Examples 



Reggeon 



r{R) 



m 



m 



m B 



Isospin + G-parity 







0/7 



V 
or 
PS 



m 



k B 



K 



K 
K* 



TT 



K 



V 



PS 



PS + V 



1.5 



k 



m B 



K 



TT 

K*{K) 
K 

7r(^) 



V 



1.6 



PS 



PS + V 



2.4 



2.0 



0.3 



Table 1: 
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Diagram 
(exchange) 



Subclasses / Examples 



Reggeon 



r{R) 



PR (GeV) 



A(*) 



PS 



2.4 



m 



k B 



else 



V 



1.6 



0.5 



0.8 



PS + V 



2.0 



B 



{ss) 



m 



A" TT 

n 



A 

K- P 

(f) n 



N 



2.8 



0.5 



1.4 



3.8 



0.3 



3.8 



0.09 



1.4 



m B 



TT p 



Table 2: 
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Ropes + 


Ropes + 


no Ropes + 




Rescattering 


no Rescattering 


no Rescattering 


p 


157.8 


199.7 


184.5 


n 


161.4 


217.6 


204.9 


A 


48.7 


35.3 


26.0 


S+ 


17.7 


12.9 


7.6 


SO 


17.8 


13.1 


8.0 


E- 


17.9 


13.3 


8.3 




5.4 


4.2 


2.1 


E~ 


5.4 


4.2 


2.0 


P 


5.6 


27.9 


11.3 


n 


5.6 


27.9 


11.4 


A 


3.8 


10.7 


2.3 


E+ 


1.5 


4.6 


0.8 


EO 


1.5 


4.6 


0.8 




1.5 


4.6 


0.8 


S° 


1.1 


2.0 


0.1 




1.2 


2.0 


0.2 


7r+ 


642.7 


692.9 


856.2 


7r° 


678.7 


724.9 


884.2 


7r~ 


680.3 


728.8 


888.9 


7^+ 


130.3 


79.0 


79.2 




75.4 


50.4 


53.0 


Ks+Ki 


200.6 


127.0 


132.1 


V 


81.7 


84.7 


86.7 



Table 3: 
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Figure 1: 




(b) 




Pi Q Ps p'j Pb P/ Ps 




(d) 




Pi p'j Pb Pf 
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Figure 2: 



(p,q) ®(l,0)= + ffi(p-l,g + l)ffi 
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Figure 3: 
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Figure 4: 
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Figure 5: 
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Figure 7: 
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{QQ)z formation 
in rope end 'plate' (b) in rope interior 
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Figure 10: 
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Figure 11: 



